A PROOF AND EXTENSION OF THE JORDAN-BROUWER SEPARA- 
TION THEOREM* 


J. W. ALEXANDER 


|. ‘The theorem on the separation of n-space by an (i | )-dimensional 
manifold} suggests the following more general problem of analysis situs. 

Given a figure C of known connectivity immersed in an n-space H, what can b 
said about the connectivity of the domain [1-C restdual to C? 

It will be shown that a certain duality exists between the topological invariants 
of Cand H-—C, and that when ( is an (7 — |)-dimensional manifold the separation 
theorem is merely one aspect of this duality. The paper also touches upon a 
number of well known related questions, -among them, the invariance of dimen 
sionality and regionality, the approachability of points of © from //-(’, the in 
variance of the topological constants of ( and //—C, and so on. 

Of course, the main difficulties in such problems as the above are of a point 
theoretic order. ‘They all yield, however, to simple pinching processes, except 


for the use of which the following treatment will be purely combinatorial. 


The earlier sections, §§ 2-S, are expository and give a rapid though essentially 


complete survey of the terminology and combinatorial machinery needed in the 
sequel. ‘The fundamental part of the discussion, with illustrations and appli 
cations, is really all contained in §$§ 9-12. 

The theory of connectivity may be approached from two different angles 
depending on whether or not the notion of sense is developed and taken into 
consideration. We have adopted the second and somewhat simpler point of 
view in this discussion in order to condense the necessary preliminaries as much 
as possible. <A treatment involving the idea of sense would be somewhat more 


complicated but would follow along much the same lines. 


* Based on a paper presented to the Society, April 29, 1916 

t First completely proved for the case » = 1 (Jordan’s theorem) by O. Veblen, these 
Transactions, vol. 6 (1905), p. 83, ef. addendum, vol. 14 (1913), p.65 (second footnote), 
Proved for a general n by L. EK. J. Brouwer, Mathematische Annalen, vol.71 (1912), 
p. 314. A proof of Jordan’s theorem along lines similar to those used in this paper was given 
by the authorinthe Annals of Mathematics, vol. 21 (1920), p. 180 
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CHAINS AND THEIR COMBINATORIAL PROPERTIES 


2. Certain advantages of symmetry are gained by setting the problem in the 


space of #-dimensional spherical geometry rather than in euclidean n-space. 
We shall therefore take as our fundamental domain the n-sphere 
+x = 1 

in the space of x + 1 real variables. The geodesics (great circles) determined 
on the n-sphere //" by its intersection with 1 — 1 linearly independent n-planes 
through the origin play the réle in this geometry of the straight lines in ordi 
nary euclidean space. A region will be said to be convex if any two of its 
points may be joined by one and only one geodesic are made up of points of the 
region. 

An n-plane through the origin subdivides the 1-sphere H” into a pair of n- 
regions bounded and separated by an (n |)-sphere. ‘The latter may in turn 
be subdivided in the same way into a pair of (7 |)-regions separated by an 
(n 2)-sphere, and so on down to a pair of 0-regions, or points. The resulting 
partition of the -sphere, consisting of two regions of every dimensionality from 
0 ton, will be called an elementary subdivision of H". It is evident that any / 
region of an elementary subdivision may be cut up by an m-plane through the 
origin into a pair of convex k-regions separated by a convex (k | )-region 
and that by repeating this process of repartitioning, the m-sphere may be cut up 
into arbitrarily small convex regions of dimensionalities 0 to 7. These regions 
will be called /-cells, where & denotes dimensionality. 

If the repartitioning is done in a perfectly random fashion, there is nothing 
to prevent the boundary of a k-cell from containing a part but not all of a cell 
of lower dimensionality. It is then possible to carry the repartitioning still 
further, beginning with the boundaries of the n-cells and working down to the 
boundaries of the 1-cells, until, finally, the boundary of each cell consists of 
complete cells only. The resulting collection of cells will be called a subdivision 
of the »-sphere H". 

A subdivision S’ will be said to be derived from a subdivision S if it can be ob- 
tained from S without the use of any other operations than the repartitioning of 
cells in the manner just described. Thus, by definition, every subdivision is 
derived from an elementary one. It should be observed that the operations 
of repartitioning may always be performed in such an order that all intermediate 
figures will themselves be subdivisions. This would be the case, for example, 
if the cells were repartitioned in such an order that no one of them was ever 
touched until the cutting up of its boundary had been completed. 

3. There is, of course, no difficulty in writing out explicitly the analytical 


expressions that determine a cell. They consist merely of the equation of H 


taken in conjunction with certain linear equalities and inequalities depending on 
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the -planes that cut out the cell. From the form of these expressions, it follows 
at once that every cell is convex. Consequently, any two z-cells are homeo 
morphic with one another (i.e., in point-for-point continuous correspondence), 
and the boundary of any 7-cell is homeomorphic with an (¢ — 1)-sphere, since 
the latter bounds each of two 7-cells on an 7-sphere. 

4. Any set of cells from the same subdivision of //" will be called a chatn 
provided the set never contains a cell /. without also containing all cells on the 
boundary of /. A chain may therefore be a very mixed agglomeration of cells. 
If, however, it consists only of z-cells and the cells of their boundaries, it will be 
called an 7-chain. 

The simplest 7-chain is one containing a single 7-cell. It will be called a 

ellular i-chain. Any 2-chain will be said to be the sum modulo 2 of the cellular 


i-chains determined by its individual 7-cells; in svmbols, 
(mod. 2) 


We shall also speak of the sum of two or more arbitrary 7-chains of the sub 
division, and it is here that the modulo 2 feature of the operation first comes into 
evidence. ‘To form the sum, we express each 7-chain in terms of its cellular 
components, as in (1), add components, and reduce coefficients modulo 2. In 
other words, the sum-chain contains the 7-cells that belong to an odd number of 
chains of the sum, but no others. 

5. An i-chain K' (7 > 0) will be said to be closed if each of its (7 | )-cells 
belongs to the boundary of an even number of its z-cells; otherwise, it will be 
called open, or bounded, and its boundary K'~' will be the (7 |)-chain deter 
mined by such (7 1)-cells as belong to the boundary of an odd number of 1- 
cells of K‘. It will be convenient to express the relation of A'~' to K* sym- 


bolically by the notation (adapted from the congruences of Poincaré) 


(2) (mod 2) 


which may be read “K’ is bounded by A’ rhe expression 


K* =0 (mod 2) 


signifies that K’ has no boundary and is therefore closed. We shall frequently 
condense the notations (1) and (2) by omitting to write mod 2. 

A 0-chain will be open or closed according as it consists of an odd or an even 
number of points. 

The boundary K‘~' of an i-cell is a simple illustration of a closed (¢ — 1)- 
chain, since every (2 — 2)-cell of k'~' belongs to precisely two (7 — 1)-cells of 
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Now, the relations 
KN! Kt = (mod 2) 


evidently imply 
(mod 2) 


or, in words, the boundary of a sum of 7-chains is the sum of the boundaries ct 
the 7-chains themselves. Thus, in particular, the sum of two or more ciosed 
?-chains is itself closed, when it does not vanish. It also follows that the bound 
ary of an open (7 +- |)-chain is a closed 7-chain, since it is the sum of the bound 
aries of the individual (7 + 1)-cells of the (¢ +. 1)-chain. 

i. We proceed to define the connectivity numbers of a chain C. Let A’ 
be an 7-chain of (, that is to say, a chain composed of cells of C. Then, by a 


second adaptation from Poincaré, we shall write 
~0 (mod 2, C) 


(A is homologous to sero, or bounds on C) provided A” is the boundary of some 
open (¢ + 1)-chain of C. The relations A, ~ 0 and Ay ~ 0 evidently imply 
Ky + Ao ~ 0; therefore, it will be legitimate to operate with homologies as 
though they were linear equations modulo 2. The expression A ~ Ky will, 
of course, be just another way of writing A’) + AL ~ 0. We denote by R*—1 
the maximum number of closed non-bounding 7-chains A5, ..., of 
that are independent with respect to homologies; that is, such that there exists 


no relation between the chains A, of the form 


R'—1 


e AS ~ 0 (mod 2, C) 


unless the coefficients ¢, are all zero, and such that every other closed 7-chain 


K' of C is related to the chains A’, by an homology 
~ Kt (mod 2, G). 


The number k’ is called the 7th connectivity number of C.* It will be seen 
later on that /’ is not only an invariant of ( but also of the set of points deter- 


mined by (. For the moment, it will be sufficient to observe that R° denotes 


the number of separate connected parts of (. It is sometimes advantageous to 


* The numbers R' are the modulo 2 analogues of the Betti numbers of Poincaré. They 
were first introduced in a paper by Professor Veblen and the author, Annals of Mathe 
matics, vol. 14 (1913), p. 163. 
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consider connectivity numbers of higher dimensionalities those that of any 
cell of C. Such numbers are automatically unity, from their definition. 

The connectivity number &’ of C which is of the same dimensionality 7 as the 
cells of highest dimensionality appearing in ( satisfies a relation that we shall 


now recall for future reference. Let there be a 7-cells 
F 


in ©, and let the symbols associated with these i-cells be regarded as variables 
free to take on either of the values 0 or |. Then, to every choice of a set of val 
ues for these variables such that at least one variable is not zero, there may be 
associated an 7-chain of ( determined by the 7-cells with symbols unity. Con 
versely, by reversing the process, to every 7-chain of ( there may be associated 


a set of values of the variables. Now, if 7 > 0, let us write a modular equation 


(mod 2) 


corresponding to each (1 1 )-cell Ay of ( ’, where the coefficient €,,, has the value 
unity or zero according as the cell Aj ' is or is not on the boundary of the cell A,. 
Then, to any set of values of the variables Aj, satisfying the simultaneous equa- 
tions /,, there will be associated a closed i-chain of (, and conversely. For, 
among the variables with non-vanishing coefficients appearing in each equation, 
there must be an even number or zero that have the value unity, which means 
that each (¢—1)-cell of C is on the boundary of an even number or none of the 
i-cells of the 7-chain determined by the solution of the equation /,. Thus, if 


p be the number of linearly independent equations /,, the maximum number of 


bes 
independent solutions must be 


which is the relation we set out to find. 


If: = 0, there are no equations /,._ For this special case, we evidently have 


We say that a closed n-chain is irreducible if its invariant R" has the value 
2, that is, if the chain is not the sum of two or more closed sub-chains. ‘The 
points of an irreducible closed n-chain form an n-dimensional manifold. 

7. Just as in the case of subdivisions, §2, a chain C’ will be said to be derived 
from a chain C if it is one of a sequence of chains beginning with C and such that 


each member of the sequence is transformable into the next one by partitioning 


a single k-cell E* into a pair of k-cells Fi and Fs separated by a (k — 1)- 


cell / . The invariants of the chain (’ are the same as those of C, for, when 
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the cell E* is cut, the only i-chains to appear in the figure that are not mere 
subdivisions of old ones are the (k — 1)-chains containing F*-" and the k- 
chains containing one but not both of F* and a The latter chains are open 
and therefore do not increase the connectivity numbers; nor do they decrease 


them by setting up new relations of bounding among old (k — 1)-chains, for 
their boundaries all contain the new cell F*~'. ‘The (k — 1)-chains containing 
F*~", whether open or closed, are transformed into old chains which no longer 
contain F*~! by the addition of the boundary of F*: therefore, they can have 


no effect on the connectivity numbers, one way or the other. Thus, we pass 
from C to C’ by a series of operations which do not alter the connectivity num- 
bers. 

[HEOoREM S”. The connectivity numbers of any subdivision S” of an n-sphere 
are all unity except the nth one which has the value 2. 

For, by §2, the subdivision S” is derivable from an elementary subdivision 
and therefore has the same connectivity numbers as the latter. But in an ele- 
mentary subdivision, every closed chain bounds a cell, with the exception of the 
#-chain determined by the two 1-cells. 

COROLLARY. Any closed (n 1)-chain K of the subdivision S” of the n- 
phere bounds exactly two open n-chains. Moreover, these two n-chains have only 
the points of K im common. 

For, since the number K is unity, there must exist an open u-chain A” such 


that 
(mod 2, C), 


and consequently a second -chain A” + S", determined by the -cells of S’ 


which do not belong to A", such that 
K" + S* = kK" (mod 2, C) 


and such that it has in common with A” only the points of A”. 
Now, if there were a third n-chain L” bounded by A”, there would be two 


independent closed n-chains 
L* + K" =0 and L" + K* + S" =0, 


and the number R” would be at least 3, contrary to the theorem. 
8. Let H” and H” be an m- and an n-sphere respectively, and let C be any 
chain of a subdivision of H”. Then, by an extension of terminology, we shall 


* This is a somewhat weakened form of the theorem on the separation of n-space by a gen- 
eralized polyhedron. For a proof of the latter theorem making use of modulo 2 equations, 
see O. Veblen, these Transactions, vol. 14 (1913), p. 65, and vol. 15 (1914), p. 506. 
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peak of any set of points C of H” in reciprocal one-one continuous correspon- 
lence with C as a chain immersed in H". For example, if C is the boundary of 
. 2-cell, C may be any simple closed curve of H". The cells of the chain C will 
ye the images in H” of the cells of C, so that C and C will both have the same 
ellular structure and, consequently, the same connectivity numbers. We shall 
requently make use of the fact that there exists a derived chain of ( made up 
That such a chain does exist follows at once from the 


f arbitrarily small cells. 
iniform continuity of the correspondence between the closed sets of points deter- 
nined by C and C, for we know that the cells of C may be redivided to any de 


‘ree of smallness. 

If the chain C does not fill up the entire space H", the residual part of H” 
will form a certain domain H” — C made up of inner points. We proceed to 
lefine the connectivity numbers of this domain. Any chain of any subdivision 
of H” will be called a chain of H” — C provided it is wholly contained in H"— C. 
Among the chains of H" — C will be set up the following homologies: (1) Each 
closed 7-chain will be said to be homologous to its derived chains; (2) each closed 
chain which bounds an open (7 + 1)-chain of H” — C will be said to be homolo- 
sous to zero. We combine homologies (1) and (2) like linear equations modulo 
2 and denote by (R’ — 1) the maximum number of linearly independent closed 
i-chains of H” — C. 
finite in this case, since we are now dealing with equations in an infinite number 


A priori, there is no reason why the number K’ should be 


It will be proved further on, however, that the numbers X’ are 
C, in spite of 


of variables. 
all finite and also pure topological invariants of the domain H” 
the fact that a metric on H” ha’ been used in defining them. ‘The number R° 
is of particular importance and evidently denotes the number of separate con- 
nected regions in H” Cc. 

Since we shall only be concerned with the relations between chains under 
homologies, it will be legitimate to do away with the distinction between a chain 
of H” — C and its derived chains. We shall therefore regard any two chains 
with a common derived chain as equivalent chains, to be denoted by the 
same symbol Kk’. A closed 7-chain will then be said to bound if it bounds in any 
of its derived forms, so that the terms bounding and homologous to zero will hence- 


‘forth be synonymous. 
ON THE DUAL CONNECTIVITIES OF C AND H” — C 


9. We now come to the body of the discussion. 

THeoreM 7°. Let C’ be a cellular i-chain ($4) immersed in the n-sphere H". 
Then, the connectivity numbers of the domain H" — C' residual to C' are all unity. 
In other words, every closed chain L* of H” — C' bounds. 

The theorem is trivial if 7 is zero, in which case, C' reduces to a point C°. 
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For there are no closed chains in H” C’ of dimensionality greater than o 
equal to n. Moreover, by the corollary to Theorem S", every closed (n — 1 
chain of H"” — C° bounds twice in H” and therefore once in H" — C°%. Ever 
closed chain of lower dimensionality bounds as often as we please in H" — C 

The general case will be handled by induction with respect to 7. We shal 


assume the validity of Theorem T’~ and first prove a lemma. 


Lemma U". Let the cellular i-chain C' be subdivided into two cellular i-chain 
A and B, respectively, meeting tn a cellular (4 1)-chain C'~' Then every k 
chain L' of H" — C' which bounds both in H" — A and H" B must also boun 
in 

If k = » — 1, the chain L* bounds exactly two open n-chains of H” meetin; 
in L*, by the corollary of § 7. The connected set of points = A +] 
must therefore lie wholly within one of these two open n-chains, since it doe 
not meet L*. Consequently, the other n-chain lies in H" — C', and L* bound 
in this region. 

Now, suppose that k<  — |, so that there exist chains in //" of dimension 
alities as high as k + 2. By hypothesis, there exist two open chains L*** and 
L" such that 

(mod 2, H" — A) 
(mod 2, H" — B) 


(4 


0 (mod 2, H” — C'~*) 


as illustrated in the figure which is purely schematic. We shall assume that 


and meet and A, respectively, otherwise the lemma would be true 


without further argument. 
Now, by Theorem 7° ", which we are assuming in the induction, there exists 
an open (k + 2)-chain VW" * ~ such that 


(mod 2. H* — 


This chain cuts the chains A and / in mutually exclusive closed sets of points 
and may therefore be broken up into cellular (k + 2)-chains so small that no 
one of them meets both A and JB, since there is a definite interval of separation 
between two non-overlapping closed sets of points. Now, let M*** be the sum 
of the cellular (k + 2)-chains of M** *that meet A, and therefore not B, and let 
L**' be the boundary of M“**. ‘Then we have 


and these combine to form the closed chain 
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ut +L ' meets neither 4 nor and therefore lies in H Hence, 


y (5) and the first relation in (4), 


pes? 4 ing (mod 2, H” 


B 


vhich establishes the lemma. 

Theorem 7° now follows at once by the ordinary pinching process. If the 
heorem were false for the chain C", it would be false for one of two cellular sub- 
hains, by the lemma, and by repeating the argument it would be possible to 
nd a sequence of sub-chains C! of C’ closing down upon a single point C° 
nd for each of which the theorem would be false. But, by Theorem T°, every 
losed k-chain L* of H” — C', and therefore of H” — C®, would bound a chain 


g of H” which did not meet C° and which therefore could not meet all of the 
chains C, converging on that point. Therefore, the theorem must be true, since 

the assumption that it is false leads to a contradiction. 


Coro.LarRy V'. A cellular i-chain immersed in an n-sphere H” cannot fill 
H", 

For let C’ be broken up into two cellular parts, A and B, as in Lemma UL", 
ind let P, and P be points of A, and B, respectively but not of the chain "li 
common to A and B. Then, by Theorem 7° ~! the 0-chain P., +P, bounds 
a l-chain in H” — C'~' which must contain a broken line of geodesics connecting 
P, with Py. But this broken line meets A and B in mutually exclusive closed 


B 
ci" A 
K+l 
L A 
\ 
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sets of points and must therefore contain points that belong to neither of thes: 
sets. Such points must be points of H” — C’. 

Corotiary W*. Let C be the sum and C’~* the intersection of two closed set 
of points A and B. Then every closed k-chain L* (k<n — 1) of H" — C whic! 
bounds a chain hy ' of H" — A and a chain Ls*' of H” — B must also boun 
in H" — C provided the chains L** and may be so chosen that Lt" + hey 
bounds in H" — 
C'~* is the null set. 


For the proof of Lemma U’ is applicable here with scarcely a change. 


Moreover, the corollary is valid evenif k = n—1 unles 


10. We are now in a position to prove the duality theorem mentioned in th 
introduction. In order to separate out the difficulties, however, let us firs 
consider an important special case which admits of a simpler proof than th 
general one. 

‘THeorem X". Let C' be ani-sphere immersed in an n-sphere H". Then the con 
nectivity numbers R* of C’ are related to the connectivity numbers R° of the resid 
ual space H"— (¢ , by the equations 


R= '=2, R=R*'=1 (s# i). 


The theorem states, in other words, that there exists but one independent 


closed non-bounding chain in H” — C’ 


This chain will be of dimensionality 
(n —2—1). It will be said to link the i-sphere 

If i = 0, the 7-sphere C" is a pair of points, so that the theorem is both trivial 
and obvious. ‘The (n—1)-chain linking C’ is any closed (x — 1)-chain of L"~' 
such that one of the points of C’ lies in each of the two open #-chains bounded by 


L"~' in H", ($7). All closed chains of lower dimensionalities bound as often as 


we please in H” — C”’. 

The case 7>0 will be solved by induction with respect to7z. Let us sub 
divide the z-sphere C’ into a pair of cellular i-chains A and B meeting in an 
(i —1)-sphere C’~*. Then by Theorem T’, every closed chain L* of H"— C’ must 
bound two open chains he and —* in H” — A and H” — B, respectively. 
Therefore, by Corollary W', L* must also bound in H” — C' unless the closed 
chain + fails to bound in H” — C’~'. But, by Theorem 
which we have a right to assume in the induction, this can only occur if 
L**! + Lt" is the(n — 7)-chain linking C’~*, in which case k = n — i — 1. 
Consequently, 


R= R= (s 1). 


On the other hand, the chain L”~‘ of H” — C'~' which does link C*~* nec 
essarily meets both A and B in mutually exclusive closed sets of points, for if 
it failed to meet A, for example, it would bound in H” — A, by Theorem 7”, 
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and hence a fortiori in H” — C'~’, contrary to hypothesis. The chain L”~* 
may thus be written as the sum of two open chains (cf. Lemma UL"), 


(6) L**= + 
lying in H” — A and H” — B, respectively, and having a common boundary 


C. 


For, otherwise, there would be an open chain L”~’, such that 


is then the required chain linking C’. 


(mod 2, H” — C’), 


and, consequently, one or the other of the closed chains L”~' + L"\~‘ or L”~* + 


L ' would have tolink C’~’, since their sum L, ‘+ L ‘would. But suppose, 


for example, that L"~'+ linked C’~'. Then, L"~' + would 
have to meet A, which would be impossible since, from their definitions, neither 
of its parts L”~‘ nor L’\"‘ could. A similar contradiction would arise if we as- 
sumed that ‘linked 

Finally, L"~‘~* is the only independent chain linking C’. For, if M"~*~' de- 
note any (7 — i — 1)-chain of H” — ci linking Cc. there is associated with 
a closed chain + ' linking C’, defined after the manner of 
the chain + in (6) associated with L"~'~'. 


Consequently, with the 
chain L”~‘~* + M”~'~*, there is associated the chain 


+ + + Mi’) 


which cannot link Therefore, by Corollary 
bounds in H” — C’, and M"~‘~* is dependent on L”~‘~*. ‘Thus, finally 
R' R*-*-! = 9. 


The chain of H” — C’ which links C’ may evidently be chosen to be irre- 
ducible, for if it consisted of several irreducible parts, one at least of these parts 
would have to link C’. 

11. This brings us to the central theorem: 

THEeorEeM Y. Let C be any chain immersed in an n-sphere H". Then, between 
the invariants R' of C and the invariants R' of H” — C there exists the following 


duality relation: 
(0<i<n-1). 


To lay the foundations for a proof by induction, let us first examine the 
trivial case where the chain C consists of 0-cells only. Obviously, a closed 
chain of H”—C of dimensionality less than n—1 bounds as often as we please in 
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H"—~C, so that R"~‘~' = R' = 1 (i>0). To determine the remaining con- 
nectivity number k"~*, let us make a subdivision of H" such that each point 
A® of C appears as an interior point of some cellular n-chain M? of the subdivision 
and such that no two of the points A; belong to the same cellular n-chain M”. 
The boundary of each cellular n-chain M! will be a closed (7 —1)-chain which 


we shall denote by 


(7) ] 2 (mod 2, H”). 


Now, every closed (n—1)-chain L"”" of H” — C is homologous to some combi- 


nation of the chains L!~°; for the chain L"”~ surely bounds in H”, 


(mod 2, H"), 


and if the bounded chain VM" contains points of (, we have merely to add to 
(S) such relations of the set (7) as correspond to the points in question to obtain 


an open n-chain free from points of C and bounded by L"””" together with some 
l 


. - n 

linear combination of the chains LL, 
Finally, there is one and only one homology between the chains L'~". For 


by §7, any linear combination of the chains L, ~ bounds exactly two n-chains in 
H", one of which must be free of points of C if the combination is to be bound in 
H" — C. Evidently, this can only occur if the linear combination includes all 


the chains 1” A Thus, if there be a points to C, we have 
= — I, 


by §6, which establishes the equality of R”~' and R’°. 

As usual, we treat by induction the case where the chain C contains at least 
one cell of dimensionality greater than zero. 

Let 2 be the chain obtained by leaving off an7z-cell of C of the highest dimen- 
sionality possible, and let A be the cellular 7-chain determined by this 7-cell. We 
shall assume that the theorem holds for / and shall prove that when A is restored, 
every change of connectivity on / is balanced by a dual change in the residual 
space so that the theorem continues to hold for C. 

Evidently, such changes of connectivity as are produced on / by the addition 
of A are caused by the appearance of new independent 7-chains containing A 
or the disappearance of independent (7 1)-chains by bounding open 7-chains 
containing A. We denote the boundary of A by C'~' and distinguish two cases 
according as 

(I) The chain C'~' does not bound on B; 

(11) The chain C'~' bounds some open i-chain A’ of B. 
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In the first case, no new closed 7-chain can be created, since such a chain would 
be of the form A + A‘ which would imply 


A' (mod 2, B), 
contrary to hypothesis. On the other hand, we have 
A GC (mod 2, C), 


so that the independent chain C’“ is lost. Noother independent chain is lost, 
for a second relation 


A+D=D'"' (mod 2, C) 
would imply 
D=C-'+pD' (mod 2, B) 


showing that D'~' was dependent on C'~'. 

The second case is treated with equal facility. Without going into details, 
we find that a single independent 7-chain A + A’ is gained and that no indepen- 
dent (¢ — 1)-chain is lost. 

To calculate the compensating changes of connectivity produced in the resid- 
ual space, we define a closed (nm — i — 1)-chain of the residual space which 
will be said to be dual to the chain A. Let an irreducible (m — /)-chain 
L"~‘ of H" — C’~' te chosen linking the boundary C’~' of A and therefore meet- 
ing A in a closed set of points. Whenever we can, we shall choose the chain 
L"~* in such a way that it contains at least one point not of A. It will then be 
possible, by the process already so frequently employed, to break the chain L”~' 


up into a pair of open chains bounded by a closed chain L”~‘~' of H" — C' and 


such that one of the open chains which we shall call M"~' contains all the points 


of intersection of L"~' with A but no point of B: 


(9) aie (mod 2, H” — B). 


‘The chain L"~‘~' will be said to be dual to A. It evidently links any 7-sphere 
- contained in C and containing A, for it has been obtained by exactly the construc- 
tion given in the proof of ‘Theorem X° for finding the chain linking such an 
i-sphere. Moreover, as we note for future reference, if « be any positive constant, 
the chain M"~' bounded by L”~'~' may evidently be so chosen that each of its 
points is within a distance ¢ of some point of intersection of M"~‘ with A. 

We derived the dual L"~'~' of A on the assumption that a chain L”~‘ linking 
fact, it is easy to prove that any chain linking 
Rather than digress to prove this, however, let us merely say that if no chain 


could be found which contained a point not of A. Now, as a matter of 


1 : 
must contain such a point. 
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/."~' contains a point not of A, then M" ‘and the boundary of M” 


is the “null” (n—7:—1)-chain. 
We now prove two lemmas. 

a) Ifaclosed chain L‘*' of H" — B does not link the boundary C'~' of A, 
there is always some chain of H" — C (that ts to say, some chain not meeting A) 
which is homologous to L**" in H"—B. If the chain L**" does link C'~', how 
ever, there is no chain of H" — C homologous to grrr. 

For if L**' links C’~' it must cut the cellular 7-chain A, otherwise it would 


bound in // A, by Theorem 7", and therefore a fortiori’ in H’ ( 


‘The same is true of all chains homologous to 1 


Qn the other hand, if L does not link ( , there exists an open chain 


such that 


AM (mod P, //" 


We may again use the figure going with Lemma l” to represent the situation 

schematically, provided we let | represent the chain L**" and imag 

ine that 1.“ \° is nearer to the eve than / and does not intersect 4, though 

may. Now M° meets A and if at all, in mutually exclusive closed 

sets of points, therefore we may repeat the argument of Lemma l” and find a 


chain / + | which does not meet 4 (nor / inthiscase). But, with the 


same notation before, 


M*"* = mod 2, HH” 


There lore, 


showing that there remains in 1 (’ a representative of the family of chains 
homologous to /. in Bb. 

Thus, in view of Theorem \, the addition of A to / cannot reduce the 
number of independent non-bounding chains of the residual space aside from 
those of dimensionality » — 7. Moreover, not more than one independent 
n — 1)-chain can be destroyed, since the sum of two chains linking (” 
cannot link 

8) Let L* be aclosged chain of H" — C which bounds in H" - 


(mod 2, H”" — B) 


n 
and, necessarily, of course,in Ho — A 


10) (mod 2, — A), 
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by Theorem T. Then, L fails toboundin — Cif and only: 


k+l 


(11) 


k+l 
+ Lp 0) (mod 2, 


links C'~" for every posible choice of the chain L'y’'. 

For, as we saw in (a), ' must cut 4 if it links Con 
sequently, L5,”' must cut A, since L\* ' does not (see figure). But if this occurs 
for every choice of L cannot bound in H’ 

On the other hand, if iy ' can be so chosen that L\” ' + 1), | does not link 
("~" then, by Corollary 11", L* must bound in H" — C. As an immediate con 
sequence of (8), no new independent non-bounding closed chains of dimension 
alities other than x ? — | can be created in the residual space when A is 
added to b. 

Let us now combine the two results obtained. 

I’) If the space //" — B contains an (n i)-chain linking no 
new independent ( 1 — 1)-chain can be created in the residual space 


lor, in place of relation (10) in Lemma (8), we may equally well write 
mod 2, // 


i — 1), and so obtain in place of (11) the expression 


But not both of the chains (11) and can link ("~~ '; therefore must bound 
in H" — C, by Corollary 1”. The net result of adding A to / is therefore to 
diminish by unity the number ” ‘ but to leave invariant the remaining con 


nectivity numbers of the residual space. 


If the space — contains no (n — 7)-chain linking ( none ot 


the connectivity numbers & can be diminished. However, in this case, a chain 
of H" — C'~' linking C’~' must meet both A and B. Consequently by a literal 
transcription of the proof of Theorem .\", we find that a single new independent 
n — 1— 1)-chain is created which is nothing more than the dual of A. 


is therefore increased by unity. It also follows that 


The number k” 
every independent non-bounding (m — i — 1)-chain of H” — C is homol 
ogous in that region to some linear combination of the duals of the 7-cells of C. 

Now, it will be observed that the changes of connectivity (1’) and (II’) of the 
residual space are exactly the ones wanted to compensate for the changes of 
connectivity (I) and (II), respectively, of the immersed figure. Furthermore, 
whether or not similarly numbered changes occur together when A is added to #4, 
the differences R' — R’~' and R"~'~! 
case and therefore remain equal to one another. Thus to complete the proof we 


° 
x both increase by unity in every 
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have only to show that the number K’ of C is equal to the number R"~‘~' of 
H" — C. ‘This we proceed to do. 

Let C’ be the chain obtained by eliminating all the cells of C of highest di- 
mensionality 7. ‘Then, to each cellular (7 — 1)-chain A’ of C’ may be found a 
dual (n — 7)-chain L”‘ of H" — C’ so close to A’ that it meets only such i-cells of 
( as contain A’ on their boundaries. ‘The latter cells it must meet, however, 
since it links their Loundaries. We may therefore break up the chain L”' 
into a set of open i-chains each containing the points of intersection of L"”' with 
one alone of the 7-cells of C and each bounded by the dual of the chain deter 
mined by that i-cell. ‘Thus, if 157 ‘~' denote the dual of A}, we shall have a 
set of homologies 


(mod 2, H" — C) 


where the coefficients €,, have precisely the same significance as in the rela- 
tions of (§ 6). Furthermore, every homology among the duals of 


chains Aj, is expressible as a linear combination of the fundamental homolo- 
gies For suppose 


(12) M*~' ~0 (mod 2, H" — C). 


Then, by combining with (9), we obtain the closed chain 
M"~' +) Mt~'=0 (mod 2, H" — C’). 
But this chain is expressible linearly in terms of the duals of the cellular (¢ — 1)- 


chains of C’so that the homology (12) is expressible in terms of the homologies 
We therefore have 


(cf. § 6), which establishes the equality of R’ and Rk" ~' 
12 
JorDAN-BRouwER THEOREM. /f M"~ be an (n — 1) dimensional manifold 
(§ 6) immersed in an n-sphere H", the residual domain H" M" consists of 
exactly two connected regions. For R® = R* = 2. 


In closing we give a few corollaries of the fundamental theorem. 


Let us distinguish between these regions by arbitrarily calling one the interior 
and the other the exterior of M"~ '. Then, 

If the manifold M"~ ' be homeomorphic with an (n |)-sphere, the connectiv- 
ity numbers of both interior and exterior are unity. For R’ = R"-'~ "=1(i>0). 
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ACCESSIBILITY THEOREM. I/n every neighborhood of every point of M" 
there is a point of M" ' which is accessible from any point P of H" — M"~' 
a broken line of geodesic arcs made up of points of H" >. 

' into cells. ‘Then, it is always pos 


For consider any subdivision of M 


sible to join the point P to a point Q on the other side of M" by a 


broken line which meets M in one (7 |)-cell of the subdivision only. Be- 


cause, if that one (2 —1)-cell were omitted, we should have k” = R® = |] and 
would become connected. But this (7—1)-cell 
may be chosen in an arbitrary neighborhood of an arbitrary point of / ; 


the space residual to M” 


which proves the theorem and also the following corollary: 

There are interior and extertor points tn every neighborhood of every potnt of 
M"~'; for example, points of the broken line PQ. 

By the same device of omitting an arbitrarily small (#7 — | )-cell, the theorem 
on the invariance of dimensionality may be proved. 

Let C* be any cellular k-chain (k < n — 1) immersed in an n-sphere H". Then, 


there are points of H" (” in every neighborhood of every point of C’; namely, 


points of a |-chain linking the boundary of the omitted (7 1)-cell. 

The theorem on the invariance of regionality is, of course, an immediate con 
sequence of the separation theorem: 

Let C" be a cellular n-chain immersed in an n-sphere H". Then no interior 


point of C” is a limit point of points of H" — C". For the boundary of C” is 
an (#7 — 1)-dimensional manifold separating //" into H’ (~ and the interior 
ac. 
Finally, we note that Theorem Y establishes the purely topological character 
of the invariants RK 


For k”~‘~' does not depend at all upon the par- 
ticular cellular structure of C, but only on the set of points determined by C; 
therefore, the same must be true of Kk’. Conversely, k’ is not affected by the 
choice of the metric on H"; therefore, neither is R"~'~'. 
PRINCETON UNIVERSITY, 
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OSCILLATION THEOREMS IN THE COMPLEX DOMAIN* 


KINAR HILLE?# 
INTRODUCTION 


The aim of the present paper is to throw some light on the question of the 
distribution in the complex domain of the zeros of functions satisfying linear 
homogeneous differential equations of the second order. The real zeros of such 
functions are well known from the works of numerous mathematicians from 
Sturm and Liouville down to living writers, but our knowledge of the complex 
zeros is very deficient. It is only in special cases that progress has been made. 


The number of complex zeros of a hypergeometric function in the case of real 


parameters has been determined by Hurwitz, Van Vleck and Schafheitlin. Hur- 


witz has also investigated Bessel functions. ‘The same field has been covered 
by Macdonald, Porter and Schafheitlin. 

Hurwitz used in his paper on Bessel functions} certain integral equalities 
analogous to those frequently used for establishing the reality of the character 
istic values in a boundary problem.§ In his thesis || the present writer used 
similar equalities for the study of the zeros of Legendre functions 

The present paper contains a systematic study of integral equalities, called 
Green's transforms, which are adjoined to linear differential equations of the 
second order. It is shown that these equalities give information concerning 
the distribution of the zeros of a function satisfying such an equation. From 
the knowledge that a particular solution of the equation in question vanishes at 
a point in the complex plane, is real on an interval, or similar information, we are 
able to assign certain regions of the plane, containing the point or the interval 


where this particular solution cannot vanish. In such a fashion we can estab 


* Presented to the Society, October 28, 1922 

+t The author wishes to express his most cordial thanks to Professor G. D. Birkhoff for 
constant encouragement and for many valuable suggestions received during the work on the 
present paper 

tMathematische An nalen, vol. 35 (1889), pp. 246-266, 

§ Cf. for instance, Bécher, Lecons sur les Méthodes de Sturm, Paris, 1917, pp. 73-76. 

Some problems concerning spherical harmonics, Arkiv fOr Matematik, Astro- 

nomi och Fysik, vol. 18 (1918), No. 17. 
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lish the existence of various sero-free regions. It is often possible by simpl 
means to make these zero-free domains cover the greater part of the plane, so 
that the zeros of the solution are distributed over a comparatively small region 

The investigation is completed by a study of the asymptotic distribution of 
the zeros by means of a singular integral equation of Volterra’s type. This 
latter method often gives information on the general disposition of the zeros and 
forms a valuable complement to the method of assigning zerv free regions 

In the greater part of the paper we treat a special differential equation, «' 
a 0, for purposes of illustration. ‘The solutions of this equation are ex 
pressible in terms of Bessel functions of the first order. 

The author has applied his methods with success to the study of several 
special differential equations suchas, for instance, Bessel’s, Legendre’s,* Mathieu's 


and Weber’s equations. ‘The results will be published in later papers. 


2. GREEN'S TRANSFORM 


2.1. Transformation formulas. We can always assume that our differential 


equation is written in self-adjoint form 


The functions G(z) and A(z) we suppose to be analytic in a region 7 of the z 


plane, where, further, A(z) does not vanish. We can replace the equation 


(2.11) by the system 


- where we have put 


* On the zeros of Legendre functions, Arkiv fOr Matematik, Astrono 


mi och Fysik, vol. 17 (1922), No. 22. 


adit’ 
— | (s) + G(s)w = 0. 
d: d: 
du’ | 
dz 
» 12) 
= 
dz 
= WwW; 
(2.121) 
= K(z) 
az 
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If we further replace the first equation in (2.12) by its conjugate, the system 


becomes 


1 = 
du, = - Wed, 
(2.122) 
dws = —G(z)udz. 


Multiply the first equation by w», the second by w, add, and integrate the 
result between the limits z; and z. Thus we get 


d- 
(2.125) ES i's + °G(z)dz = 0. 


Here we have assumed 2, and 22 to be points in 7 and also the path of integration 
to be entirely in 7. This expression (2.13), which is of fundamental importance 
for our coming work, we shall call the Green's transform of the differential equa- 
tion on account of its similarity to a well known formula of Green for real vari- 
ables. 


If we write 


= dK, aa 1dK> = dK, 


G(z)dz = +idPr. = 


we have 
(2.15) “dK + 2dr = (), 


or, if the real and the imaginary parts be taken, 


(2.152) Ym ES + + = (). 


2.2. Change of variables. Let us introduce a new independent variable in 
system (2.12) by putting 


(2.21) dz = f(Z) dZ, 


dz 

— 
(2.14 K(2) 


w 
w 
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whereby the system is transformed into 


du wdZ, 
(2,29) k(Z) 

dws 2(Z)wdZ, 
with 

K(z 
k(Z) = 4 g(Z) =G(z)f(Z). 
{(Z) 


Z Z 
(2.23) ES —_ 2 dk + f = 0, 
Z Z 


where 
- = dk, 
k(Z) 


(2.231) 


2(Z)dZ = dy. 


Thus we infer that the Green’s transform is invariant under a transformation of 
the independent variable. 
By special choice of Z we can obtain simpler forms of the transform. Our 
first choice will be 
(2.24) Z = K(:z). 
If we put 
J(Z) = G(z)K(s), 


we obtain 
Z 2: Z 
2.241) ES 2dZ + 21(Z)dZ = (0), 
Zi Zi 


to which formula corresponds a differential equation 


(2.242) + = 0. 
This special type of differential equation is the one with which we shall be con- 
cerned mostly in this paper. 
We can of course also choose 
(2.25) Z = T(z),* 


* This function has of course nothing to do with the classical '-function. 


= 
The Green’s transform goes over into 
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3: 
which yields 
2 


with 


Finally we can get a more symmetric form by putting 


? 26) dZ ad 
\ 


which leads to 


d7 
S(7) 


with 
S(Z) vG(s)A (2). 


The differential equation takes the following form 


252) 4 
7, dZ 


This transformation we shall use in the last section of the paper for the study of 
the solutions in the neighborhood of an irregular singular point. 

Up to this point we have changed only the independent variable. We can of 
course also make a simultaneous transformation of both independent and de 
pendent variables 

f(Z), 


Starting from an equation of the form (2.242), for instance, we can find in 


finitely many transformations of type (2.27), which furthermore preserve the 


form of the differential equation. For each set of such variables Il’, Z we have 
a Green's transform of type (2.241) of our equation. 


The simplest possible of all such transformations is 


W, 
o= aZ + b, 


(2.28) 


TUT | / * / dZ (), 
(/{ ) (,( VAL ( 
corresponding to the differential equation 
> ad | du | 
(4) t % (). 
dZ dZ 
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which carries /(z) into a*?/(aZ + }). We shall return to this transformation 
$3.8. 
It is evident that w; and we play the same roéle in the system (2.12 From 
this remark it follows that we have also formulas like 


, d 


This is of course nothing but the conjugate of the expression in formula (2.15 
We shall find later on that in all theorems concerning distribution of zeros, « 
and enter symmetrically. 

2.5. Twonets of curves. Wecan evidently simplify the formulas (2.151) and 
(2.152) by choosing the path of integration so that one of the following relations 
is fulfilled, namely 

dK, dK» (); 
dP, (); (0. 


These paths constitute four families of curves which together form two in 
dependent nets of orthogonal trajectories. 
In the case of an equation of type (2.242) the K-net is simply the net of lines 


parallel to the axes. The same lines serve as a P-net of an equation of type 


On account of the skew symmetry between G(z) and A(z) we need consider 
only one of the two nets. Our results will hold, mutatis mutandis, for the 
other net. We choose the P-net, as the K-net is trivial in the most import 
ant of all special cases, namely that in which A(z) l. 

Through every point a of the region 7 in the s-plane where G(s) is regular 
and furthermore G(a) # 0, passes one and only one curve of each of the families 
Ir, and M, which curves we denote by ya and Pea, respectively. These two 
curves are, of course, orthogonal to each other, with slopes equal to g)(a) gela 
and g(a) gi(a), respectively, if we put 


+ 


Hence the F\-family has horizontal tangents along the curves g;(z) 0 and 


vertical ones along go(z) = 0. For the P.-curves the state of affairs 1s reversed. 
We shall now proceed to consider the behavior of the I'-curves in the neigh- 
borhood of an exceptional point of G(s) where this function cither vanishes or 
becomes infinite.* We restrict ourselves to the case where G(z) is analytic in 
7 except for poles. 
* Cf. for the theory of curves of this kind F. Lucas, Géométrie des polynémes, Journal de 


l’ Ecole Polytechnique, vol. 28 (1879), and F. Klein, Ueber Riemann's Theorie der 
Algebraischen Funktionen und threr Integrale, Leipzig, 1882. 


(2.22). 
(2.052) 
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First, let z = a bea k-fold zero of G(z). Then 


k+1 
G(z)dz = a, (s + 
a k + l 


i¢ 


0 
a=re anda, we have 


p +1 aos [« 
| 


Hence (k + 1) branches of Iya pass through z = a and the same number of 
branches of Ta intersect them there in such a manner that the tangents at a 
of the one curve make equal angles with each other and bisect the angles formed 
by the tangents of the other curve. 


If s=a is a pole of order k > 1 of G(s), we have 


With the same notation as above we obtain 


l 

l 


where y +76 = C. Let us assume for the sake of simplicity that no logarith- 
mic term appears. Then we can conclude that every curve of the two families 
which passes in the neighborhood of z = a actually passes through this point. 


The I\-curves are tangent to the lines 


arg (0 a) = [e+ +1) (k—1) (v=0,1,... ,k —2). 


Replacing 2v + 1 by 2» we get the tangents to the P.-curves. 


In the case where z = a is a simple pole of G(z) we have 


Gea: = C +a, log (z —a) + 


2 \ 
G(z)dz ¢ = 
J ( 
/ G(2)d: =( — 
k—1 (s—a) 
| / | =y- ; 
om | / =6+- 
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Lf =y+a logr—B,6+ 
Jim 6+ 6, logr+a,0+ 


where y + 62 = C, and a; + Byt = a. 
If neither a nor ;, is zero all curves of the net in the vicinity of s = a admit 


and 


this point as an asymptotic point which they approach in the manner of a log- 
arithmic spiral. When 6, = 0 all curves of I, sufficiently near to a are ovals 
around this point which are approximately circular for great values of y. In 
the first approximation, the I',-family behaves like a pencil of lines through a. 
If a, = 0 the réles of the two families are interchanged. 

When carrying over the results of the discussion of the descriptive properties 
of the I'-net to the K-net, we have to observe that the rdles of zeros and of 
poles are interchanged; a zero of A(z) being a pole or logarithmic point of K(z) 
and so on. 


2.4. G(z) a polynomial. Suppose G(z) isa polynomial of degree 1 
G(2) = Ao(z — (2 — ae)” ... (2 — a,)”™ =n). 


A non-specialized curve belonging either to I, or TF; has no double points at 
all in the projective plane. Such a curve does not pass through any of the 
points a), a2,'°*, a,. There is, however, one curve of each family which passes 
through the point a, and has a (vy, + 1)-tuple point there. Thus there are at 
most yw singular curves in each family. 

The intersections with the line at infinity in the projective plane are all dis- 
tinct for one and the same curve but the same for all curves of the family. The 
asymptotes are all real and distinct; they all intersect in a fixed point, the center 
of gravity of the roots + + + v,a,)/n, where they furthermore 
form equal angles with each other. If arg Ay = ¢» the asymptotic directions 
of a Ia are given by 


arg 2, = 


_ 
n+1 2 


The asymptotes of the [:-curves pass through the same point and bisect the 
angles between the asymptotes of the I'\|-family. We note that the asymptotes 
of the curves g:(z) = 0 and ge(z) = 0 also intersect in this point. 

A I-curve can never begin or end in a finite point, nor can it be closed. This 
implies that every point in the plane can be reached along a suitably chosen 


|| 
n+ 1 
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path from infinity without crossing the curve in question. This is a conse 
quence of the fact that P\(z) and P.(2), being harmonic functions throughout 


the finite plane, can not have maxima or minima. 


3. DETERMINATION OF ZERO-FREE REGIONS 
The real axis. We shall 
4.1, ame real axis. We shall apply the results of the preceding section to the 


problem of locating the zeros of a given solution of a linear differential equation 


of the second order. We start with the real zeros and restrict ourselves to the 


case when A (2) From formulas (2.151) and (2.152) we obtain 


Now take an arbitrary solution w of the equation (2.242) ora pair of solutions 


of the corresponding system. Can this solution w, or, more generally, 


can the product w; @» vanish at two points 1; and v» of a certain interval (a, > 


of the real axis? Suppose g keeps a constant sign in (a, } Then by (3.12) 


Jey Can vanish once at most in the interval in question, and similarly if ¢, 


\ < 


yn (a,? We have assumed of course that G' /(z) is analytic in the inter 


val. Thus 
THEOREM 33.1 l of the real axts throughout which is analytr 
| keeps aconstant sign (perhaps vantshing 


wid furthermore 
is an arbitrar\ 


of & T 


Thus if wis to oscillate on the axis of reals, supposed free from singular points, 


we must have g;( > 0 and g(z) must change sign or vanish identically. 


In case there are singular points of the real axis we may be able to get addi 
Suppose for instance that ais 


and A» such that 


tional information concerning the zeros 
a regular singular point of the equation with exponents A, 
MA»). Then /(s) has a pole of at most the second order 


+ re > 
at a. Tf we assume WA) > | 2 the expressions in (3.11) and (3.12) will 


remain finite when we let +, or x converge towards a, provided w, is the solution 


corresponding to the exponent A,; and 0 when x» —* a. Conse 


quently this particular solution will not vanish in an interval (a,)), a being the 


singular point, where ¢g)(z) and go(z) are subject to the conditions stated in the 
theorem above. This remark, of course, also applies to more general types of 
singular points. 

It is obvious that the results of this section can be generalized to equations of 


the type (2.11). 
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5.2. Dynamical interpretation. Theorem ‘}.1 admits of direct dynamical in 
terpretation 
If we put u u + iv in equation (2.242), assume 


separate reals and imaginaries, we obtain the following system 


These are the equations for the motion of a particle of unit mass, subjected 
to a force whose components on the radius vector and on a line perpendicular 
to the same are g,; vue? 4 and goVu- 4 , respectively. 

Theorem 3.1 now gives a sufficient condition that a particle starting from the 
origin in the w,7-plane at the time x, with a certain velocity shall continue to 
move away from the origin during a certain space of time. The condition 

v) S 0 evidently means that the radial component of the force is always 
directed from the origin whereas the invariability of sen means that the 
rotating component has the same sense during the time interval in question 
These conditions are evidently sufficient to prevent the particle from turning 


~ 


towards the origin. 


3.3. Linear paths. We shall now consider an arbitrary linear segment 


}in the complex plane and an arbitrary solution w(s) of equation (2.242 


under the assumption that /(2) is analytic along the segment chosen. Put 


From (2.241) we infer 


If we put 


cos 24 
we obtain 
dw 
dr 


1y22 
+ re". 
ar + ¢ hd) () 
ad? . 
d 
0S 26 vols) + sin ()(s, 
dr + / P(s, 1); 


360 EINAR HILLE [June 


These formulas yield the following 


THEOREM 3.31. There is at most one zero of the product w dw/dz, where w(z) is 
an arbitrary solution of w" + J](z)w = 0, on a segment (2, 22) in the complex plane 
where ](z) is analytic, provided either 

(i) P(z, 0) <= 0 along the segement, or 

(ii) O(z, 6) keeps a constant sign there. 

If in addition to (i) f\(O) 2 O, or in addition to (ii) f2(0) has opposite sign to 
that of Q(z, 0) all along the segment, then there is no zero at all of the product w dw/dz 
on (2, Ze). 

Let us now extend this theorem to a region 7 in the complex plane where 
J/(s) is analytic. ‘Take a pencil of parallel lines /, 


Zo? 
— re 


Each line is characterized by the variable point z. Let 7 be simply connected 
and of such a shape that every line of the pencil cuts the boundary in two points 
at most. Then there will be two lines of /,, each of which meets the boundary 
in two coincident points, a and 8, respectively. ‘These two points a and 8 divide 
the boundary in two parts, one of which we shall designate by C and which shall 
be the locus of the points z. Then we have 


THEOREM 3.32. Jn a region 7, defined as above, there is at most one zero of the 
product w dw/dz on any line lz, if either 

(i) P(z, 0) = 0 throughout T, or 

(ii) O(c, 9) 4 O in T. 

If in addition to (i) Rle“w(z)dw/dz| = 0 along C, or in addition to (ii) 
3 m[e"w(z)dw dz|= 0 or has opposite sign to that of Q(z, 6) in T all along C, then 
there 1s no zero of the product in T. 

The first part of the theorem does not lay any serious restrictions on the zeros 
of w, and the second part seems rather unmanageable. We shall however give 
some special cases of great interest. Suppose C is a segment of the real axis, 


8 
@ = m/2 and w(z) takes on only real valuesonC. Then Rle"wdw/dz]= Oand 


we obtain 


THEOREM 3.321. If the equation (2.242) has a solution w(z) which is real on a 
segment (a, b) of the real axis; if further T is a region symmetric with reference to 
(a,b) such that every vertical line which cuts the region cuts its boundary twice and 
meets (a,b) in an interior point; and if finally R{J(z)| = O throughout T, then 
neither w(z) nor dw/dz can have any complex zeros in T. 

If in the statement of the preceding theorem we replace the real axis by the 
axis of imaginaries, vertical lines by horizontal and finally suppose #{/(z)] < 0 
in 7 the theorem still holds, the exceptional zeros now being purely imaginary. 


1922] OSCILLATION THEOREMS 361 


As a simple application we may choose the equation w” + w = 0 with w = 
sin z. Theorem 3.321 shows that there are no complex zeros of sin z or cos z. 


We can of course also derive a special theorem from formula (3.33), namely 


THEOREM 3.322. If the equation (2.242) has a solution w(z), real on a segment 
(a,b) of the real axis where w dw/dz keeps a constant sign, namely the same as that 


of 3m[J(z)] in a region T above the real axts, whose base ts (a,b) and of which the 


boundary above the real axis 1s cut in one point at most by any vertical line; then 
neither w nor dw/dzcan havea complex zero in T or in the symmetric regionT below 
the real axis. 

A similar theorem can be stated for the imaginary axis. 

As a simple illustration take the equation 


dz?" g 


This equation will be used throughout the remainder of the paper for the purpose 
of giving concrete applications of our theorems. ‘Thus we need some knowledge 
of the nature of itssolutions. % isan irregular singular point whereas 0 is regu- 
lar singular with exponents0 and 1. There is one solution, = vz _/\(27vz), 
where /;(a#) denotes the ordinary Bessel function of order 1, which is an entire 
function of z and vanishes at the origin; every other solution is of the form 
lc: + ce log z|E(z) + ceH(z) where H(z) is an entire function. (z) is real for 
all real values of z and has infinitely many real negative zeros; every other solu- 
tion can be real on a half-axis at most and if it is real on the negative half of the 
real axis it must oscillate there. A solution which is real for x > 0 can have 
at most one positive zero or extremum. In general w® will increase beyond all 
limit when x —> + *%. In fact we have 


w'(x)vx 


= +1, 


(3.351) lim 


+H w(x) 


save for one exceptional integral which yields —1 in the limit.* 

As # (—1/z) > 0 when R(z) = x < 0 we infer by theorem 3.321 that no 
integral of the equation which is real for real negative values can have any com- 
plex zeros in the left half-plane. Further {m(— 1/z)>Owhenx >0. ‘Take 
a solution which is real for real positive values of z and let x» be the positive zero 
of w(z) dw/dzif there is any; otherwise positive but arbitrary. Then by theorem 
3.322 we infer that such a solution can not vanish in the half-plane W(z) = x». 

* Cf. A. Wiman, Uber die reellen Lésungen der linearen Differentialgleichungen zweiter Ord- 
nung, Arkiv f6r Matematik, Astronomi och Fysik, vol. 12 (1917), No. 
14. 


d*w l 
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lor the exceptional integral we may in fact state a little more, namely that it 
can not vanish in the right half-plane including the imaginary axis. 


lor these exceptional integrals we may give a general statement. Let 


(3) 
where / (2) is real and positive for real values of s > «+; and analytic in a region 
’) including the real axis for (zs) > x,. Concerning /) we make the further 


assumption that any horizontal line which cuts /) shall cut its boundary in one 
finite point only, and along the part of the line which lies in /) one of the two 


relations 
MNLF (s) | > () or 


shall hold Let II be a solution of 


such that O when s — ~ in/) along a parallel to the real axis. Such 
an exceptional integral will under very general assumptions on /'(z) tend to 
zero so rapidly that the integrals in formulas (3.33) and (3.34) converge when the 
upper limit tends to x. This will for instance be the case when F(z) is a poly- 
nomial. We conclude that such an exceptional integral can not admit a finite zero 
in D nor can its derivative vanish there. 

The results of this section can be considerably extended in various directions. 
‘lo every theorem we have stated for equations of type (2.242) corresponds a 
theorem for equations of general self-adjoint form. The reader can easily work 
out these theorems for himself. 

4.4. The star. In the preceding section we have considered parallel lines. 
Let us now take a pencil of lines through a finite point z = a where /(z) is regular 
and does not vanish. Put 

L(s) = (s — a)? /(s) = P(s) + 1Q(s), 


and draw the curves P(z) = 0 and O(z) = 0. ‘These curves intersect at s = a 


where each of them has a double point, the tangents of which are given by 


cos 26 g(a) — sin 26 g(a) = 0, 


cos 26 g(a) + sin26 g(a) = 0, 


respectively. Let us follow a line in general position from z = a. We start 
with definite signs of P(z) and U(2), and travel along the line until O(z) changes 


sign. If perchance P(z) has been positive or has changed its sign, we break off 


() 
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the line at that point. In case, however, P(z) has been negative all along the 
segment, we continue to follow the line until finally even P(z) changes sign, at 
which point we stop. We always stop at a point where /(2) is singular no 
matter what the signs may be. In the same manner we treat every line of the 
pencil and determine on each ray a last point p,. The assemblage of all seg- 
ments (a, p,) form together a configuration which we shall call the star belonging 
toa.* ‘The boundary of the star is made up of the points p, and parts of cer 
tain critical rays which are to be found among the tangents that can be drawn 
from the point z = a to the curves P(z QO and O(z) = 0. 

The points p, are either points of these curves or singular points of /(z) 

The appearance of the zero-free star is somewhat different from that of a 
star of convergence in as much as it in general is composed of two different re- 


gions which only touch each other at z = a. (See the figure below.) In the 


Pic. 1. 


neighborhood of this point the boundary of the star is made up of that branch 
of the curve O(z) = 0 which approaches z = a through a region where P(z) > 0, 
and of the tangent to that branch at z = a. ‘The fundamental property of the 
star is given by the following 

THEOREM 3.4. = ais a sero of w dw dz where w(z) is a solution of equation 


(2.242), then this product does not vanish throughout the star belonging to z = a, 
including the regular paints of the boundary. 
This theorem is an immediate consequence of Theorem 3.31.4 


* It is, of course, the similarity in form and generation to the star of convergence, introduced 
in analysis by Prof. Mittag-Leffer, which suggests this name. The chief property of our star 
is indicated by the more explicit name zero-free star of w(z) with respect toc = a 

+ For further developments of the properties of the zero-free star ef. the paper quoted in 
section 3.9 below, also Convex distribution of the zeros of Sturm-Liouville functions, B wu 1- 
letin of the American Mathematical Society, vol. 28 (1922), 


pp. 261-65 and A correction, ibid., vol. 28 (1922), p. 462. 


C=0 P=0 
\ ZS 
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Let us use this theorem for proving that the entire function /(z), considered 
in the preceding section, as well as its first derivative, has only real and negative 
zeros. E(Q) = O, but the origin is a singular point; we easily find however 
that the integrals in formulas (3.33) and (3.34) converge when 2, is allowed 
to approach 0. Hence we can form the star for the point z = 0 and apply 
Theorem 3.4. We find L(z) = —z and the star consists of the whole plane, 
except for the negative part of the real axis, which proves our assertion. What 
we have been proving is, of course, only that /,(#) has only real zeros, which is 
a well known fact in the theory of Bessel’s functions. 

3.5. The primary domain. [et us separate reals and imaginaries in formula 
(2.241). We obtain 


(3.51) w dx + w|*(gidx — gody) = 0, 
(3.52) dy + jw|*(gidy + godx) = 0. 
dz Ja dz 


Assume = is a zero of w dw, dzlying in a domain 7 where g;(z) amd _ go(z) 


are different from zero and /(z) is regular. In order to have a definite case to 
deal with we suppose 
gi(s) > 0; < O 

in 7. Draw parallels to the axes through 2) and fix the attention on the second 
of the resulting quadrants of 7. ‘The points of this region are of two different 
kinds; a point @ may be joined with z) by a curve which never increases with 
increasing values of x = N(z), or a point may be such that there is no such 
curve. The points of the first kind form a certain simply connected region D2 
where the subscript 2 denotes the quadrant. Similarly we find a certain region 
PD, in the fourth quadrant, the points of which all can be reached from 2 by a 
continuous, never-increasing curve. ‘These two regions together form a region 
D2(29) which we agree to call the primary domain of 2. In view of (3.51) and 
(3.52) we obtain 


THEOREM 3.5. If there is a solution w(z) of equation (2.242) such that 


w(z)dw/dz vanishes at a point z in a region T where J(z) is analytic and 
ML (z)] 4 0, Jm[/(s)] ¥ O, then neither w(z) nor dw/dz can vanish in the 
primary domain of 2. 

The primary domain is a region Dy; if Sm[J(z)] > O, but a region Dey when 
Sm[/(z)] < 0. 

The primary domain gives of course only a first rough approximation of the 
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zero-free region in the neighborhood of a given zero, but even that is useful for 
orientation. 

We have assumed that 2 is not a point of a curve g;(z) = Oor g(z) = 0. If 
2) however lies on such a curve, we can still find a similar zero-free domain. 
Suppose, for example, that a segment (a,b) of the real axis forms a part of go(z) = 
0 and that a < z < 6; then we can form a primary domain which is of type 
Dy, if go(z) > 0 above the segment and of type D2; in the contrary case 

The boundary of a D,,(zo) is made up of parts of the curves gi(z) = 0 and 
g2(z) = 0 and of lines parallel to the axes. It is as a rule possible to cross these 
boundaries and continue the zero-free region beyond the limits of the primary 
domain. The formulas (3.5) indicate clearly the way of doing this; we shall 
take up this question in a more general form in the next section. 

3.6. The standard domain. Let us take our differential system in the general 


form 
du l 
dz K(z) 
(3.61) diz = — 
dz 


(3.631) dv = dV, + = G(z)dz, 
dz 
(3.632) dK = dK, + 1dK2 = 


We assume G(z) and A(z) to be single-valued analytic functions of z, existing 
in the whole plane. Then the two nets fF and K are uniquely determined; 
one and the same curve may however correspond to different values of the con- 
stant of integration, depending upon which branch of the, in general, many- 
valued functions ['(z) or K(z) we are considering. 

A solution of (3.61) is as a rule a many-valued function which undergoes a 
linear transformation when z describes a closed contour in the plane surrounding 
singular points of the system. These are the zeros a, dz, ..., @,, ... of A(z), 
the singular points by, be, ..., b,, ... of G(z) and K(z) and, in general, the point 


with the corresponding Green’s transform, split up into real and imaginary parts 

21 | 21 
(3.622) | + |w»| 2d Ke + 120, = 0), 

21 21 
where as usual 
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at infinity. In order to have to deal with single-valued functions only, we can 
either join the singular points of the system by properly chosen cuts which ¢ is 
not allowed to cross, or we can construct a Riemann surface with branch points 
at the points (a,,), (/,,) and * on which every solution of (3.61) is single-valued. 
We use the latter scheme and imagine the P-net and the K-net represented on 
the Riemann surface. 

Now take a point a on a particular leaf of the surface where we know that 
W"(s) wy(s)uy(c) = 0. In order to show that another point > on the surface 
cannot be a zero of Iz) we try to join a and } by a suitable path of integration. 
Formulas (3.62) suggest the use of the two nets for this purpose. We define 
a path curve as a standard path if it fulfils the following conditions : 

i) /t does not pass through any singular point; 

ii) /t vs composed of a finite number of arcs of curves belonging to the two bas. 
nets and K; 


iti) -V/l along the path one and the same of the followitne four characteristic in 


equalities ts satisfied, namely 
ar, < 0, adr, = 0, dV. = 0, dV. 0, 
(33 1) 
dK, = 0: dK, < 0; dK. = 0; dK» < 0. 


With the definition we have chosen the tangent of the path is continuous ex 
cept ata finite number of points where ares of different base-curves mect. We 
make the additional agreement that in a small neighborhood of such a point 
the two meeting ares are replaced by an are which joins continuously and has 


a continuous tangent. This can always be done in such a manner that the 


characteristic inequality of the path is preserved. 
According as the path in question is characterized by the first, second, third 
or fourth set of inequalities, we call the path a standard path of the first, second, 


third or fourth kind, respectively 
SK,, SK,, SK. and SK, 


Take a point a on the surface where G(s) and A(*) are analytic and A(a) # 0 


and construct all standard paths which start at a. In this fashion we obtain a 


certain continuum of curves spread out on the surface. The set of points which 
belong to at least one of these paths forms a region on the Riemann surface 


which we call the standard domain /)(%,) of We include the boundary points 


of 1)(s)) in the standard domain except the point 2» itself and eventual singular 


points of the differential equation. The fundamental property of /)(29) is given 


by 
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‘THEOREM 3.6. /f cy ts a sero of W(s) = wi(z)wels) then Ws) cannot van 
ish in the standard domain of 

We note that the primary domain serves as a kernel of the standard domain 
in case A(z) = 1. 


In our typical example we have 


K(z) - log s. 


The K-net consists of the lines parallel to the axes; the P\-family is made up of 
the circles with center at the origin and the P.-family of the straight lines through 
the same point. The solutions of (3.55) are single-valued on the Riemann 


surface of log 


plane 
Q 
| z= 
| 
t 
Fic. 2 
Take a solution such that where x) < 0. A simple discussion 


shows that this function must be real for real negative values. If we construct 
the standard domain /)(x») we obtain the whole plane (except the origin) co. 

ered once and, in addition, two overlapping regions Wt(s) 2 lvul, >0 and 
Wis) 2 xo, Yin(s)< 0, attached to the lower (upper) side of the positive real 
axis. Using the method which is explained in the subsequent section, we could 
even prove a little more, namely: every solution which is real on the negative 
real axis has only real negative zeros in the sector — (32,2) << arg(—2)< (jm 2 


Now take a solution for which IV’ (x%») = O with x» >0. Such a solution must 


he real for real positive values of s. The standard domain of x» consists of 


ty 
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(1) the half-plane (zs) 2 x» in one leaf of the surface; 

(2) the region < xy, R(z) Vand | 3m(z)| xo, Owhere the negative 
part of the real axis is counted as a double line; this all in the same leaf as 

(3) the quadrant (zs) < 0, Ym(s) > 0 in another leaf, hanging on region 
(2) along the lower side of the real axis; 

(4) the symmetric quadrant in a third leaf, hanging on region (2) along the 
upper side of the double line. 

If the zero is not real we get standard domains which are asymmetric. ‘The 


figure gives the typical shape of )(a + 1b) when a< 0, b > 0. 


3.7. The standard domain of acurve. In the preceding sections we have 
formed various zero-free domains of a point. We can easily extend this con- 
ception to an are of a curve. 

Let ( be an are of a continuous curve on the surface, having a continuous 
tangent except at a finite number of points. We assume further that G(s) is 
analytic along ( and A(z) analytic and different from zero. Let a solution of 
(3.35) be given along (. By given we understand that there is some means at 
our disposal by which we can ascertain the signs of the real and of the imaginary 
parts of w (dw dz) (zs). 

‘Take a point 2%) on the curve (. From 2» emanates a set of standard paths, 
belonging to all four families. Of these families we can as a rule only use two 


as paths of integration. Put 
(3.71) w(z) goA(s) = +2 V(2). 
ae 


If, for example, ('(2o) > 0 we can use SK, but not SK,. ‘Thus we draw 
from every point on ( the standard paths of the kinds which are indicated by 


the table below. 


| negative | zero | positive 
I 
(3.72) | SK, | SK,,SK, | SK, 
Vis) | | SKs,SK, | SK, 


| 


These paths together form a continuum of curves whose points make up the 
standard domain of C with respect to w(z), which we denote by DCw. The 
points on the boundary of (Cw are included in the standard domain by defi- 
nition except the points on ( and the singular points of the differential equa- 


tion. In view of formulas (3.62) we have 


THEOREM 3.7. Jf a solution, w(z), of (2.11) be given, in the sense defined above, 


along a curve C on the Riemann surface of w(z), which curve does not pass through 
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any of the singular points of the differential equation, then there is no zero of wdw/ dz 
in the standard domain of C with respect to w(z). 

We have assumed that C does not pass through any of the singular points of 
the differential equation. This condition can, of course, be abandoned in 
special cases with great advantage. If, for instance, z = a is a regular singular 
point with exponents A, and dA» where ‘(A,;) > 1/2 say, then the expressions in 
formulas (3.62) remain finite when z, is allowed to approach a, if w;(2) is the solu 
tion corresponding to Ay. Further w)(z)w.(s) —» 0 when a. Conse- 
quently the standard domain of such a point exists and has the same properties 
as the standard domain of an ordinary regular point. 

In case s = a is an irregular singular point we are often able to find certain 
sectors with vertex at a and corresponding solutions which tend toward 0 when 
> — > a within the sector and approach zero so rapidly that formulas (43.62) are 
applicable with a as the lower limit of integration. In this case we can not con- 
struct a full standard domain around a, as we are allowed to use only such stand 
ard paths as, in a small neighborhood of a, lie on the sectors mentioned above. 
Such sectors can in general be found when, for instance, the irregular singular 
point is of finite rank. 

3.8. Change of variables. In section 2.2 we have shown that the Green's 
transform is invariant under a transformation of the type 


where /(z) is, of course, an analytic function of <. 

Take a point % in the s-plane and form the corresponding standard domain 
(so). Suppose F(z) is analytic in D and F’(z) ¥ 0 there. Then Z = F(z) 
will map /)(z9) conformally on a region d(Z») of the Z-plane. We assert that 
d(Z») is the standard domain of 7) in the Z-plane. In fact, every standard path, 
SK, starting from 2%» in the z-plane, is transformed into a curve, sx, that starts 
at Z,) = F(z») in the Z-plane. In view of the formulas in section 2.2, 


(3.81) d T(z) dy(Z); dK(sz) = dk(Z). 


Hence the base nets in the two planes correspond to each other. From this it 
follows that sx is a standard path in the Z-plane, starting at 7» and, furthermore, 
sk is of the same kind as SK. Thus every point of d(Z») belongs to at least one 
standard path sx; consequently d(Z») is a standard domain as asserted. We 
can express this fact also by saying that the standard domain of a point is an abso- 
lute covariant under conformal transformations of the plane. 

All the different forms of the Green’s transform, given in section 2.2, are 


consequently equivalent for our purpose; one choice of independent variable 


f(Z) or Z = F(z) 
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may however lead to simpler expressions for k(Z) and g(Z) or have some other 


advantage over another choice of variable. In the case of our typical equation 


the transformation in formula (2.25) offers some advantages. By putting 


Z=-—logz or 


d Z dw 


very solution is now single-valued, in fact, an entire function of 7. This 


the equation becomes 


fact more than compensates the difficulties arising from the slightly more com 
plicated form of the x-net. 

In order to get any new facts concerning the zero-free regions we must change 
both variables simultaneously. We shall not develop any generalities along 
this line. A simple example may be sufficient to show what use one can have of 
a linear transformation. 

In section 3.6 we determined the standard domain corresponding to a solution 
of our typical equation with a positive zero or extremum. Tet us indicate the 


hypothetical) complex zeros of the principal branch of the solution in question 


by = , and Sas where as usual z, stands for the conju 
gate of s,. Puts, = x, + 7). Our result showed only that from a certain 
value of n, x, << x» and y, > x9. We shall prove in section 4.7 that y, = 
Oly x,). Thus it is desirable to obtain an upper limit for y,, for all values of 1. 


This can be done as follows. Put 


eZ. 
Equation (3.35) becomes 
d*x 
Thus we have 
= — e” log Z. 


. Putting Z = ke we get 


yi(Z) = sin 0@ — cos 6 log k, 


= — sin R — cos 
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which shows that the y-net is made up of logarithmic spirals. The point a 


is transformed into xe ; the corresponding y,- and y2-curves are 


,ran Wo + 
R = 

, cot + 0 
R = xo ’ 


respectively. The first curve has a vertical tangent at 


It is easy to see that this point, as well as the points on the same vector for which 
R > Ro, belongs to a zero-free region. In order to determine the locus of [ty,¢, 


we introduce polar coérdinates in the old s-plane by putting > = re*. Then 


¢g = 26,r = FR and we obtain 


Thus there is no zero of w(z) dw dz beyond this curve. Hence 
Vu < ave / log 


or 
arg 2, 


n 


3.9. The domain of influence.* ‘Ihe zero-free domains considered in the 
preceding articles have all been obtained by specialization of the path of in 
tegration in formula (2.15). We can arrive at a more general type of zero-free 
region that embraces the previous ones as special cases in the following way. 


The formula 


3.91) ES ES £ kK 4 / = 0 


contains four terms. ‘The first term will evidently not vanish if the sum of the 
other three terms is different from zero. This sum will certainly not be zero if 
the argument of each of the terms in the sum is known to lie in one and the same 
angle of opening less than 7. ‘The argument of the second term in (3.91) de 


pends essentially upon the initial conditions. We have 


(3.92) arg (wwe) = arg (dw/dz)/w + arg N\(:) 


* This article was added at the revision of the paper in November, 1921 


” 
d= 06, Ro = 
r=. ,¢ tatio * 
oe = 
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where (dw/dz)/w at z = 2, is a known number \ which may be 0 or «. Put 


arg + arg A(z;) + = u, 


(3.933) — argdK = x,, 
= 
Given an angle 6 such that 
(3.931) 
we can always choose the path of integration, starting from z = 2, so that 


6<K«K,<0+7, 
0+ 7, 


at all points z on the path C, say. Then it is obvious that the arguments of the 
third and the fourth terms in (3.91) fulfill the same condition as x, and y, do. 
Consequently the sum of the three latter terms in (3.91) is different from zero 
for this particular choice of the path of integration and 

(3 94) O 

for any point 2(# 2) on C. We have of course tacitly assumed that the curve 
( does not pass through any of the singular points of the differential equation, 
an assumption that will be made throughout the present section. 

To a given angle @ corresponds a set of paths of this nature which we call the 
lines of influence with regard to the point 2, and the solution w(z|z, 4) where 2, and 
\ in an obvious manner denote the initial conditions. Giving 6 all possible 
values, consistent with (3.931), we obtain all lines of influence belonging toz = 2, 
and the solution in question. ‘These lines we imagine to be traced on the Rie- 
mann surface mentioned above in section 3.6. The set of regular points on the 
surface which can be reached by lines of influence from the point z = 2, forms 
a zero-free region, the domain of influence of 2, and w(z|z, 4) which we denote by 
DI(2|X). 

It is evident that we obtain the most extensive domain for a given point 2, 
when = Oat = 2, i.e., when Ais 0 or ~w, there being no restriction (3.931) 
in that case on 6; and the regions we obtain when \ # 0 or « are interior parts 
of DI(z, 0). Further we can verify that the standard domain and the star of the 
point z, are both contained in DI(z\0). 

In a similar manner we can define the domain of influence with regard to a 


curve and a solution.* 


* For further developments of these ideas see the author’s paper, On the zeros of Sturm- 
Liouville functions, Arkiv f6r Matematik, Astronomi och Fysik, vol. 


16, No. (1922), 17. 
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4. THE ASYMPTOTIC DISTRIBUTION OF THE ZEROS* 


4.1. Reduction to the normal form. In order to investigate the distribution 
of the zeros of a solution of a given differential equation in the neighborhood of 
an irregular singular point we use the transformation in formula (2.26). We 


have 


d dw 
(4.11) K(2) + G(z)w = 0. 


By putting 


(4.12) Z(z Z0) = f dz, 
20 


we obtain 


d dw 
(4.13) a S(Z) i + S(Z)w = 0, 
where 
(4.131) S(Z) = G(z) K(z). 


As a parallel form of (4.13) we also use 


d*w dw 
(4.1 t) dZ? + F(Z) a? + w= (), 
where 
(4.141) F(Z) S(Z)]. 


Let us further introduce a new independent variable 
(4.15) W = VS(Z)w. 


The differential equation then becomes 


where 
Z) = 1 dF 1 F(Z)) 
(4.161) = 3 [F(Z) 


* This section was rewritten at the revision of the paper. For the results compare with the 
memoir of P. Boutroux, Recherches sur les transcendents de M. Painlevé et l'étude asymptotique 
des équations différentielles du second ordre, Annales de I’ Ecole Normale Su- 
périeure, (3) vol. 30 (1913), pp. 255-375, and (3), vol. 31 (1914), pp. 99-159, especially 
$§2-6. Cf. further R. Garnier: Sur les singularités irrégulieres des équations différentielles 
linéaires, Journal de Liouville, (8) vol. 2 (1919), pp. 99-200, especially §28. 
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The function 7(2 %) depends upon an arbitrary parameter 2 and, the value 
of 2) being fixed, Z is in general infinitely-many-valued. We assume that we 
can find a determination of Z such that the function ¢(Z) is single-valued and 
analytic in a region A of the Z-plane, extending to infinity and having the fol 
lowing properties: 


CONDITION A 


1) Ais simply-connected and smooth; 
2) Every line parallel to the real axis cuts the boundary F of A either (1) 


in a line-segment, or (11) in a point, or (ili) not at all; 


3) A lies entirely in a sector 


A region A that fulfils a condition A is said to be of type A. In the case 
when Fis actually cut by every parallel of the real axis we say the region is of 
type Aa, otherwise of type Ab. From the condition A it follows that A contains 


a strip 3) given by some equality like 
1.17) M(Z7) =A>R; Be = 


Furthermore we assume that ¢(Z) in A satisfies 


CONDITION B. 


M 


Zits 


where vy and M are positive constants. 

1.2. An integral equation. ‘The method of successive approximations shows 
that the general solution of (4.15) is bounded in any strip Ay of finite width. 
Consequently the expression 


(4.21) f(Z) = WZ) + sin (7 — Z) @ (7) W(7) dl 
Zz 


has a definite meaning and represents an analytic function in Ay. Here I1)(Z) 


is a solution of 


I1(Z) is a solution of (4.15) and the path of integration is a straight line parallel 


to the axis of reals. Moreover we find that 


P"(Z) + f(Z) = (Z)N (4). 


te 


r+6* arg Z 6: sia k 
}.22) W,” 
(4.63) 
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Consequently if /(7) is a solution of the integral equation. 


+ f sin (/ Z)\o(1) fl) dT, 
then /(Z) is also a solution of (4.15). The equation (4.24) is a singular integral 
equation of Volterra’s type which we can solve in the following way. Putting 
1.2.9) (Z,T) sin (/ — Z) @(1), 


we define the iterated kernals 


T) = KAZ. at, 

KAZ. 7) = U)K,-(U, T) dv, 


and further 


(4.252) (Z) - K,(Z, 1) WT) dT. 
Z 
Then 
(4.253) W(Z) = + E,(Z) 
n l 


formally satisfies (4.24). If Ao is a strip of finite width and 7 and Z are points 
in Ay, we have 
Sh, 


sin (7 Sk, 
(4.26) 


and consequently 
(4.261) E,(Z)|< » 


Hence the series in (4.253) is absolutely and uniformly convergent in Ay 
and represents an analytic function of 7 in that region, which function is further- 
more a solution of the integral equation (4.24) and consequently also of the dif- 


ferential equation (4.15). 


379 
| 


376 EINAR HILLE [June 


We know in advance that (7) is bounded in Ay and we find easily that 


(4.27) < 2K, 
when 
(4.271) \Z| > ‘] in Ap. 
Vv 


Using formula (4.24) we can improve this result and obtain a second approxi- 
mation, namely 


(4.272) —Wy(Z)| < 


when Z is subject to the condition in formula (4.271). 
This shows that I!’(Z) approaches 11)(Z) indefinitely in A», or, using the 
terminology of Boutroux,* the function I1(7) is asymptotic to 11)(Z) in Ay. 
Similarly we find 
273 — Wy'(Z)| < 


when Z satisfies (4.271); this shows that I1’’(Z) is asymptotic to I1y’(Z) in Ap. 

1.3. Oscillatory solutions and truncated ones. By formula (4.253) we obtain 
to each solution |! (7) of (4.22) a solution I’(Z) of the given differential equation 
(4.15) and I!'(Z) is asymptotic to I1)(7). Taking a pair of linearly independent 
solutions I1y,(7) and I(7) of (4.22) we obtain a pair of linearly independent 


solutions and of (4.15). In particular we choose 
(4.31) Wo(Z) = WelZ) =e 


The corresponding solutions and are asymptotic to and 
respectively, and consequently admit of no zeros in Ay beyond a certain limit. 
‘The same is true of their derivatives. Even in the full region A these solutions 
do not admit of any zeros outside of a sufficiently large circle. We notice from 


formula (4.24) that 


= 1 + 
(4.311) a 
WAZ) = 1 


where 6,(Z)| < C, a certain constant, when Z| > K in A, which proves the asser- 


tion. Using a term coined by Boutroux7 we call these integrals truncated in A. 


* Loc. cit., pp. 270-273. 


+ Loc. cit., pp. 261-263 
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If the region A is of type Aa (vide supra) these twointegrals are the only 
truncated onesin A. In fact, every solution of (4.15) can be represented in the 


form 


(4.32) W(Z) =aWi(Z) + aW(Z), ce. 


which is asymptotic to I/o(Z) =cie"” + ce “in A. This region being of type Aa 
we can always choose the strip A, b; < Ym(Z) < Be, in such a fashion that it 
contains the zeros of W)(Z), a, = a+ mum, from a certain 1 = No on. As 
cco # O we can assume W)(Z) = sin (7 — a) without any loss of generality. 
Let us further denote the parts of A which lie above and below Ap respectively 
by A, and A-;. In these regions 


(4.330) W(Z)= sin (Z— a) + 2? Zin Ao; 
6,(Z) 

(4.331) et “W(Z)=e+t™ sin (Z — a) + = , Zin Ay; 

(4.332) sin (Z — a) + 
Here 
(4.333) (A=—1,0, + 1) 
when 

1 

(4.334) \Z| = [2M 


in A,, where A, denotes the maximum of |e" sin (Z7—a)| in A). 

Let us mark the points a, in Ao and surround each of them by a small circle 


Z—a,,| =e. On the circle T',, we have lsin (Z — a)| > ~e. If we choose 


1 
(4.34) \Z| > 


< =e. Hence the first term on the right hand side in (4.330) 


dominates the second one when Z is a point on a circle T’,, in the part of A that 
lies outside of the circle Ip defined by (4.34). We can assume without con- 


siderable loss of generality that Ip cuts the boundary T of A in two points only. 


Then Ip divides A in two parts, one interior A— and one exterior A+. We 
denote by Ax what is left of A+ after punching out the interior of the circles I,,. 

Then, W(Z) has one zero and only one tn each of the circles T,, in A+ and no 
zeros at allin Ax. ‘The first part of the statement is a consequence of a well 
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known theorem due to Rouché.* ‘The latter part of the assertion depends upon 
the fact that sin (7 — a) is the majorating term in the expression for W(Z) in 
Ax, provided ¢€is small enough. In the same way we find that the zeros of W’(Z) 
in A+ all lie in the circles T’’,,, — a’,| = e where a’ + namely one in 


n 


each circle. 

‘Thus in the case of a region A of type Aa the zeros of II'(Z) are approximated 
by the zeros of a certain sine-function, © sin (Z — a), to which IV/(Z) is asymp 
totic, and the zeros of 11’’(Z) are approximated by the zeros of C cos (Z — a). 
Let the zeros of ’(Z) in A be denoted by A,, and those of 11’’(Z) by A’, where 


n 


the notation is so chosen that 
lim (A, — a@,) Q and lim (A’,, — a’,) = 0. 


Then we say that the points A,, form a siring of zeros in A and similarly the points 
A’,.— There are only two linearly independent solutions which do not have 
a string of zeros in A, namely the truncated solutions defined above. 

In the case of a region A of type Ab, however, we have infinitely many in 
dependent solutions which are truncated in A, namely all for which the zeros of 
the asymptote-function |1)(7) lie outside of A. On the other hand, given a 
point set a, = @ + um which lies in A from a certaim 7 on, we can always find 
a solution I1(Z) of (4.15) which is asymptotic to sin (7 — a) in A and the 
zeros of which form a string in A, approximated by (a@,,). ‘This solution is given 
by our previous formulas (4.25) with 11)(Z) = sin (7 — a). 

In order to emphasize the dependence of the solution I1'(7) on the parameter 
a we denote it by I1'(7, a). The zeros of 1(Z, a) in A are functions of a@ which 
justifies the writing of A,(q@) instead of simply A 


Let us see how A,(a@) 


varies with a. We put @ = o + 77and let o increase from oy to oy) + 7, keeping 


r constant Then A,(a@) describes a certain curve from A,(ao) to 
A, (ao + 7) A,,4,(@), Which shows the relationship between the functions 
A, (a). The path followed by A,(@) is almost a straight line from a point close 
to a, a + um toa point close to a, ,,. If ¢ continues to increase, A,,(a@) 


will describe a certain curve Y(7)) that joins the different zeros in the string, 
belonging to 11(Z,a9). We call this curve a zero-curve (of solutions) of the differ 
ential equation (4.15). The curve is evidently asymptotic to the line Qin(Z) 

ry and is uniquely characterized by its asymptote. ‘Through every point in 
A+ there passes one and only one zero-curve. 


1.4. Varied conditions. In the preceding articles we have assumed ¢ (7) to 
* Rouché, Mémoire sur la série de Lagrange, Journal de I’ Ecole Pol ytech 


nique, cahier 39 (1862), pp. 217-218 
t Boutroux and Garnier both use the term /igne de z(ros for what we denote by string of 


zeros 
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be analytic in an infinite region in the right half-plane. The case when ¢(Z) 
is analytic and fulfils a condition B in a region A in the left half-plane of type A 
(with suitable change of condition A(3)) can be treated in a similar manner. 
We find that the general solution of (4.15) is asymptotic to some sine-function 
W.(Z) in A and there are solutions the zeros of which form a string asymptotic 
to a given point sect a_, = @ — urin A. 

A special case of large importance is when (7) is analytic and fulfils a con 
dition B in an upper half-plane T, defined by Qin(Z) > 42. In this case every 
solution is asymptotic to some sine-function I1),(7) in T,, the extreme right 
part of fT, and asymptotic to some other sine-function I!y4(7) in T_, the extreme 
left part of T. Thus the solution may be truncated in T, or in T_ or in both of 
these regions. If the solution in question is oscillatory in T, then the zeros form 
one and only one string which is approximated by a point set. (q@,). And if the 
solution is oscillatory in T. then the zeros form a single string that approaches 
a point set (8_,,) indefinitely. 

On the other hand given a sine-function I1o,(Z) there exists a solution 
W’, (Z) of (4.15) that is asymptotic to Il, in T,, and there is another solution 
IW (Z) that is asymptotic to Wy, in T.. But ifr =Qm(a@) is a very large 
number then there exists a solution I1’,(7) that is asymptotic to sin (Z a) 
throughout T. The zeros of Il’,(7) are approximated by @ + nm and conse- 
quently the two strings in T, and in T_ join into one doubly unlimited string. 
There are evidently no zeros above this string and only a finite number below it 
in T. 

We can of course treat a lower half-plane in a similar manner. 

1.5. The distribution in the z-plane. Let us now try topass back to the 
plane. We will be chiefly concerned with the case of a region Ain the right 


half-plane. The function 
(4.51) Z 
' / \ KX (z) 


establishes a one-to-one and conformal correspondence between the points of 
the interior of the region A in the Z-plane and the interior points of a certain 
region /) of the z-plane. ‘his latter region is of course simply connected but may 
overlap itself infinitely often. In the case of overlapping we consider )) as part 


~ 


of a Riemann surface D over the z-plane on which we proceed to study the in- 


tegrals w(z). “These solutions are of course single-valued and analytic in the in 


terior of D (on D ) which does not contain any of the singular points of the 
differential equation (4.11). ‘There is however at least one singular point on the 
boundary of ) corresponding toZ = %. 


By this conformal mapping, the zero-curves = (7) of the differential equation 


(4.15) are transformed into a set of curves S(t) in D. The points a, are carried 
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over into points a, in ); similarly the points A, into points A,. The circles 
r’,, go over into closed contours C, around a,. Finally the circular are I’, in 
A is carried over into a curve C, in ) which separates this region into two parts 
D, and D_ corresponding to 4, and A_ respectively. The points of D, that 
remain after leaving out the interior of the closed contours C,, we denote by Dx. 

‘To the solution W(Z,a) of (4.15) corresponds a solution w(z,a) of (4.11) by the 


relation 


1 
(4.52) w(z,a) = [S(Z)]~2W(Z, a). 
If a, lies in A from a certain 7 = N» on, then the corresponding points a,, will 


finally lie in D. The solution w(z,a) will have precisely one zero within 
each of the curves C,, in ),, namely at A,, and there will be no zeros at all in 
Dx. If D corresponds to a region A of type Aa there are two and only two lin- 
early independent solutions which are truncated in /), namely those corre- 
sponding to 11",(Z) and I1’2(Z), and every other solution oscillates infinitely often 
in , When the corresponding region A is of type Ab, however, there are 
infinitely many solutions which are truncated in ). 

The zero-curves of the equation (4.11) as far as the region /) is concerned are 
precisely the curves S mentioned above. ‘These curves are asymptotic to the 
curves © satisfying the differential equation 


(4.53) vm dz =(). 
K(s) 
The integral curves of this equation we shall call the asymptotic zero-curves of 
(4.11). 
Through every point ain ) (on D ) passes one curve S and only one. Follow 
this curve S(a) in the direction of increasing values of 3i(7) and mark the 


points a), dx, ..., a,,...defined by the equalities 


(4.54) dz (n=1, 2,...), 


the integral taken along (a). 
By our previous investigation we know that there is a definite solution w(z,a) 
of (4.11) having the points a,, as limiting points of its zeros in the sense that the 


actual zeros A,, can be so numbered that lim (A, — a,) = 0. Let us draw two 
circles 'p, and I'p, in the Z-plane with radii R, > Kk; > R. Suppose for the 
sake of simplicity that these circles cut I’, the boundary of A, each in two points 
only. The maps in the z-plane of these circular arcs are two curves Cpr, and Cp, 
in D which together with the boundary of D enclose a curvilinear quadrilateral 
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{[D] on D. This region is cut in two points only by © (a), 2; on Cx, and z on 


~* 


Cr, say. Then the number of zeros of w(z, a) in [/)] is given by 


4.5; ji = az - 
(4.55) N[D] ff 


where the integral is taken along © (a) and —1 S$ @ S + 1. 
Let us now consider the zeros of w’(z,a). We have 


dz |dZ 2 <4) 4) dz" 


The second and the third factors do not concern us. We note that I1'(Z) is 
bounded in A» and assume that 


M 


: (u > 0).T 
VAL 


(4.561) F(Z)| < 


Further, the function dll’ dZ vanishes once in each of the circles ['’,, and not 
at all in A‘. Using the theorem of Rouché, quoted above, we conclude that the 
function dlV/dZ — 1% F(Z)\V(Z) vanishes once and only once in each of these 
circles and is different from zero in every finite point of A+, provided the quan- 
tity « which enters in formula (4.34) is small enough. Hence dll /dZ vanishes once 
and only once in each of the corresponding closed contours C;, in D+ and not 
at all in 

Just a few words about regions A and T considered in section 4.4. A region 
A of type A in the Z-plane corresponds to a region L in the z-plane. The zeros 
form a simple string just as in the case of a region 7. A region T corresponds 
to a region U in the z-plane; this region consists of three parts V+, U— and 
U’° corresponding to T+, T— and the intermediate portion of T respectively. 
The zeros form simple strings in U+ and U- (one string in each at most); in 
certain cases these strings may join to one doubly unlimited string. 
* Compare the result in the real case given by Wiman, loc. cit., p. 7. Professor Wiman 
has informed me that he has used the more or less classical transforination in formula (4.12 
in his lectures at the University of Upsala for studying the properties of solutions of linear 
differential equations of the second order in the complex plane. The general character of the 
distribution of the zeros in the neighbourhood of an irregular singular point has been known 


to him for some years. 

+ The function F(Z) satisfies a Riccati equation (4.161). Using this equation it can be 
proved that F(Z) will satisfy (4.561) in A with » = 1 or v where v is the exponent entering in 
condition B, provided v = 1. When 0 < » < 1 the general solution of (4.161) may admit 
strings of poles in A. Leaving out properly chosen regions around these poles, condition 


(4.561) is still satisfied 
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1.4. Polynomial coefficients.* Let us assume that A(z) and G(z) are poly 


nomials 
(4.61) K(s) G(s) = + 
Infinity is an irregular singular point if g — k 2 — 1, whichis the only inter 


esting case in this connection. 
In this case Z(s\<0) is in general an abelian integral of the third kind. In the 
neighborhood of infinity we have, putting g — k + 2 = m, 


m log z ) 


(-£.62) = ae ~ get {’ 


where the logarithmic term occurs only when m is an even number. Conversely 
we find 


2 
1.1 


2V 2% 


the series being convergent in any sector in the neighborhood of infinity. Fur- 


ther we have 


G(z)N"(s) + (2) 
(4.63) 


3 
\G(z) [AX (s) |2 
Expanding in ascending powers of > we obtain 

j x 


2 | 


2V 
Hence, substituting the expansion for z from formula (4.651) 


ble series of > same 
1.632) F(Z) double series of the same 
m ZN type as in formula (4.621) { 


The expansion for @(7) is of similar form but starts with a factor 
[(3k — g + 4) 4Am?](1/Z?). 


Thus @(Z) satisfies a condition B with v = | in any region outside of a 
sufficiently large circle in which the argument of 7 is bounded. Let us take 


as our region A the part of the 7Z-plane in which |Z, > Kk and W(Z) 2 0. If 


* Cf. J. Horn, Uber die irreguliren Integrate der linearen Differentialgleichungen 2weiter 


Ordnung, Acta Mathematica, vol. 28 (1909), pp. 171-201, especially p. 198 
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R is sufficiently large this region is conformally mapped by 2(7|Z)) on a smooth 


. (yg) . 
sectorial region /)”’ in the z-plane in which 


l 
(4.64) args < + 
m m 
where 6 is a small number, #% = arg g) and yp = 0, |, ., m — 1 depending 


upon which determination of ¥/Z we use in formula (4.621). The determina 
tion of the logarithm does not affect the form of the region appreciably. 

The asymptotic zero-curves in the z-plane are given by 
m 1 
0 


(4.65) r2 sin > (md + = const. = re’) 


showing the leading term. They form m pencils in the neighborhood of in- 
finity which point they approach in the directions 


(4.651) = [ | (u=0, 1, ...,m—1). 
m 


In general w( :la) is not single-valued in the neighborhood of infinity. Let 
D be the part of the Riemann surface of log z which lies outside of a large circle. 
Then the solution is single-valued on D. ‘The zeros of w(z| a) are seen to form 
m strings asymptotic to the directions @“’ in each leaf of D. In particular if the 


solution is uniform in the neighborhood of infinity we find m strings in all. 


If N(r) denotes the number of zeros in a string within the circle z = r and 
we put 
9 m 
(4.66) = 
ma 


then N(r)/n(r) —> | whenr—~> ~. 

Let us take up the question of truncated solutions. By the results of section 
1.3 we know that there are exactly two solutions which are truncated in one of 
the directions 0”. Hence the total number of truncated solutions is at most 
equal to 2m. This number, however, actually reduces to # as we shall see. 

Let us take a modified region A bounded by a large circular are |Z R, 
z —52argZ2— + 6(6 a small fixed positive quantity) and the tangents 
at the end-points of this arc, extended to infinity in the left half-plane. 
This is a region of type A in which (7) satisfies a condition B. The two solu 
tions, and 1W2(Z), truncated in A are asymptotic to and '”, respec 
tively, in this region. Moreover they keep their asymptotic form in a wider 
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region, being asymptotic to e'” when —x < arg Z < 2x and W.(Z) being 
asymptotic to e~ '“ when —2x < arg Z < 7, as can be proved by considering a 
region A symmetric to A. Hence W’\(Z) —~> 0 in the upper half of ‘A and 
W.(Z) —»> Oin the lower half of A provided IS m(Z)|—> *. ‘These solutions 
are furthermore uniquely determined by this property, as the solution of the 
integral equation (4.24) is unique. 

‘The region A corresponds to m different regions Din the z-plane in which 


(4.67) (gp = 0,1, ...,m — 1). 


These regions )’ of course have common parts; in fact, D“’ hasa part U™’, 
< arg: 0"" _ ¢ in common with * and so on. By formula 
(4.52) we conclude that there is one solution truncated in )“’which tends 
towards zero in U™ anda second solution truncated in D“ * ? which approaches 
zero in the same region l””’. There is, however, only one solution which tends 
towards 0 in LU“. ‘Thus we find the number of truncated solutions reduces to 
m; these solutions are uniquely characterized by their asymptotic properties: 
the solution w,(z) which approaches zero when z —> © in U™’ is truncated in 


+1) 


the adjacent directions @“’ and 6“ It is worth while observing that the 
truncated solutions preserve the same asymptotic representation in three ad- 
jacent regions U'"’, namely U“-), U™ and U“*” in the case of w,(z). This 
number # of truncated solutions is actually reached as is shown by considering 
the differential equations w” + z‘w = 0. 

Let us take the case when A(z) 1. ‘Then the general solution of the differ- 


ential equation which is an entire function of z can be represented in the form 
(4.68) w,(z) = w,(z) — 


where w,(z) and w(z) are linearly independent solutions and A is a complex 
parameter. When is this solution truncated? The answer is: \ has to be one 
of the asymptotic values of the meromorphic function X(z) = w;(z)/we(z), or, in 
other words, \ has to be a transcendental singular point of the inverse function 
=(\). A quantity a is said to be an asymptotic value of an entire or meromorphic 
function f(z) if there is a path L (a Jordan curve) tending to infinity along which 
f(z) tends towards a. ‘The condition is necessary because, if w,(z) is a truncated 
solution, then we can find paths L along which w,(z)/w(z) —> 0. Hence 
w(z)/we(z) —»> along L; thus \ is an asymptotic value of A(z). The con- 
dition is also sufficient. If \ is an asymptotic value, w,(z)/we(z) —> 0 along 


the corresponding path L. But 


(4.69) w (2) /we(z) =W,(Z) /W2(Z) ~ W(Z)/Wo(Z), 
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with obvious notation. This shows that the path can be taken outside of the 
sectors6” — 6 < arg z < 6“ + 6 and, moreover, that the limit is zero only if 
W is either e“ ore But then W(Z) is truncated in A and consequently 
w,(z) is truncated in a corresponding region 2)’. 

This result shows further that every function of the meromorphic family 
[cw (2) + cowe(z) ]/[cswi(z) + cywe(z)] has at most m asymptotic values. 

4.7. Concluding remarks. The discussion in this section has given us in- 
formation on the distribution of the zeros in the neighborhood of an irregular 
singular point, in as much as we have means by which to determine the number 
of strings of zeros and their directions, and furthermore to find the restrictions 
imposed on the absolute value of the difference of two consecutive zeros in a 
string. This is often enough to give a fairly good picture of the general char- 
acter of the distribution. ‘The results in the third section, on the other hand, 
give information about zero-free regions around a given point or curve, and thus 
about the restrictions on the argument of the difference between two zeros. 
Hence the two methods are complementary. 

The reader may find it instructive to work out for himself the distribution 
problems for the differential equation 


(4.71) w’ — = 0, 


by combining the two methods. We find Z = 27 Vz and the transformed normal 


form is 
(4.72) Ww’ + a)! 


The asymptotic zero-curves of (4.71) are parabolas with focus at the origin and 
the negative real axis as axis. ‘The distribution of the zeros in the vicinity 
of the asymptotic parabola is very regular. 
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ON CERTAIN RELATIONS BETWEEN THE PROJECTIVE THEORY OF 
SURFACES AND THE PROJECTIVE THEORY OF CONGRUENCES* 


BY 
FRANK EDWIN WOOD 
1. INTRODUCTION 


Wilczynski has shown} that the system of partial differential equations 


Vy = M2, 2y = Ny, 
(1) Vun = av + bz + cy, + dz, 
= a'y + + 
where 
oy oy 
vy, = ow’ Vy = dv’ = au? etc. 
and where wm, 7, a, ...,d’ are functions of u and v, will be completely integrable 


if the following integrability conditions are satisfied: 
c=f,,@’=f,, b= —d,—df,, a’ =—c', —c'f,, mn —c'd = fur, 
Mun + dy + + — fym, = ma + db’, 

2n,n + mn, = a, + fymn + a’d, 

mn + 2mn, = b’, + fymn + be’. 


In such a case, (1) will have precisely four linearly independent solutions 
K) _(K) , ' 1 1 1) 
(K =1, 2, 3, 4). Let ands be inter 


preted as the homogeneous coérdinates of the two points P, and P,. As u 
and v vary, these points will describe two surfaces S, and S, and the line P,P, 
will generate a congruence whose focal surface consists of the two surfaces 
S, and 


* Presented to the Society, under a different title, March 26, 1921. 

t E. J. Wilezynski, Sur la théorie générale des congruences, Mémoires Publiés 
par laClasse des Sciences de l’Académie Royalde Belgique, 
Collection en4,° ser. 2, vol. 3 (1911 This paper will hereafter be cited as the Brussels 


Paper 
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The most general transformation which will leave the set (1) in the same 
form is 
(3) y = Muy, 2 = u(v)z, = alu), v = 


where \, u, a, 8, are arbitrary functions of the variables indicated. The coeffi- 


cients m, n, a, ..., d’ resulting from this transformation are 
m= n= n, = ¢ = 
AB, May hay, MB, 


to which may also be added 
1" Kef 
(. 
= 


where K is an arbitrary constant. 

There are various covariant configurations connected with the congruence. 
Some of these are connected with the focal surface and present themselves at 
first referred to a local coérdinate system connected with that surface. Others 
appear directly in their relation to the local coérdinate system of the congruence. 
In order that we may study the relations between these covariant configura- 
tions, their equations must be related to a common tetrahedron of reference. 
The fundamental covariants of the congruence form a tetrahedron, and the 
relations which exist between this tetrahedron and the local tetrahedra of ref- 
erence for S, and S,, respectively, must therefore be obtained. The develop- 
ment of these relations constitutes the first part of this article. 

There are two asymptotic curves which pass through a point P, on S,, and 
the osculating linear complexes of these two curves have a linear congruence in 
common. ‘The two directrices of this congruence are called* (d;') the directrix 
of P, of the first kind; this directrix lies in the tangent plane of S, at P, but does 
not contain P,, and (d,”) the directrix of P, of the second kind; this directrix 
passes through P,, but does not lie in the tangent plane. 

The curves cut out on S, by the parametric curves of the congruence form a 
conjugate system; the osculating planes at P, of these two curves intersect in 


* E. J. Wilczynski, Projective differential geometry of curved surfaces (Second Memoir), 
these Transactions, vol. 9 (1908), p. 95. 


] - Au Nun - l Au Quu 
(4 1 ae + ) b hae (1 “4) (. 
on c’ b’ + a’ qd’ = (a’-2 ), 
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a line which passes through P, and is called* (x) the axis of the point P, with 
respect to the conjugate system. 

The dual of the axis of P, is called{ (7:) the ray of P, with respect to the 
conjugate system. 

There are four lines d's, d"s, x2, re similarly connected with P, on S, making 
eight lines in all. In the second part of this paper, the equations of these eight 
lines referred to the tetrahedron of the congruence are obtained. 

In the remaining parts of this article the possible coincidences of two or more 
of these lines are considered in general, and congruences which are characterized 
by certain pairs of coincidences are obtained in a canonical form. ‘The existence 
of congruences which have no one of the eight associated lines indeterminate 
and which possess any one of twelve possible coincidences is established. In 
particular a study is made of the congruences which possess one of the following 
coincidence pairs: (1) the directrix of the first kind and the ray of P,, coincide, 
and the directrix of the second kind of P, coincides with the ray of P,, or (2) 
the axis of P, coincides with the ray of P,, and the ray of P, coincides with the 
axis of P,, or (3) the directrix of the first kind of P, coincides with the directrix 
of the second kind of P,, and the directrix of the second kind of P, coincides 
with the directrix of the first kind of P,. This class of congruence has a special 
interest since it can be connected with the theory of functions of a complex 
variable. 

The writer wishes to thank Professor Wilczynski for his kindly interest and 


valuable suggestions. 


2 RELATIONS BETWEEN THE TETRAHEDRA OF REFERENCE OF THE FOCAL SHEET 


AND THE TETRAHEDRON OF THE CONGRUENCE 


The differential equations which characterize S, are} 


Vuu — Ven = + CVu + 


mmm 
(5) m 
= mny + — 
m 


Make the transformation 


= (u,v), 1 = v), 


* E. J. Wilczynski, The general theory of congruences, these Transactions, vol. 16 
(1915), p. 314. 

t Loc. cit., pp. 317-18. 

t E. J. Wilczynski, Brussels Paper, p. 44, equation 79. 
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where . d 
m 

The new parametric curves “; = const. and 7; = const. will be the asymptotic 


curves of S, and the transform of (5) will be in the intermediate form.* 
Let F'(u, v) be a general function; then 
OF OF OF oy 
Ou Ov; Ou 
_ OF _ dF dy | dF dy 


whence 
1 
F,, = (F, + 
(7) | 
F, (—F,-+ 6F,) 
Using (6) and (7) one obtains 
Vuu = unr + (dp x, + Ou; 
(8) + Ver» 


Vuv = = + (bg + — (Pv + Vor» 
+ ain + + ur + Vivo 


If these values (8) be substituted in (5), we find the differential equations of S, 
referred to its asymptotic lines. These equations are 


+ uy + ay = 0, 


(9) , ; ; 
Youn + 20,Vu + +¢,y =0, 
where 
= 1 — 26— + 366, + 6, I, 
46°¢, m m m 
Wy b Ms My 
2h) = 533 — ic —-— — & — — 26 — — 66, — 5, }, 
m m m 
1 
= ws (a + 26mn), 
“(10 
_ Pv 6 b 62 26 Mu 66 6 
2b; — (452% +5 +26 —" + 388, — 
(a — 26mn). 


*E. J. Wilczynski, Projective differential geometry of curved surfaces (Second Memoir), 
these Transactions, vol. 9 (1908). 
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Let us introduce the notations* 


l , Mw bx ' ld, 3m, 
— f — log 6m?) = ——- — 

l 


1 /b mt 6: 1 oO + lop Sn? 3d My 
t\d m 6 Ov 4\"" 2d + 2m 
(11) 
Fy B, ( ly 
Gy = By, ( 
0 


then equations (10) become 


( — — (a+ 26mn), = — 3 (a — 


The fundamental semicovariants of S, aret 


(13) l ; 
Fa by Vu, + aay (a1) + J; 


and the fundamental relative covariants of the congruence aret 


My n 


From (1), (7), (13) and (14) the following relations are obtained: 


(15) 50, m 


l My 4 m l 


Py + Rip +Sioc , 


Il 


* Cf. E. J. Wilezynski, Brussels Paper, p. 20, equation 17. 

t E. J. Wilczynski, Projective differential geometry of curved surfaces, these Transac- 
tions, vol. 8 (1907), p. 248. 

t E. J. Wilczynski, Brussels Paper, p. 23, equation 23. 


l My 6 ; | 
6 om 6 
(12) 2b = “SG, 2%, == }, 
v,\v 6 
l | 
| 
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where 


Muy p 
mM 


+ + + 


(16 p Nn, 
= m —60,+2- +2 
( n 
= — 6 — 2h, , 
OP AN l 


40.~,5; =2m. 


Solving (15) for 1, 2, p and o gives 


(17) 
p= — oO Aig, — by + + 69,21 — OY 
om 
Rip + Sy, 
where 
20 ang, + by, , m 
om 
9 
< d >’ 
O; (- VY; — 0¢ Ri) 
m 
(1S) 
R, = — v (- Qi + 
WwW m 
») 
— 


From equations (11) and (12) we find 


ag, + bi, — 2C; +2) 


(19) 


| 
| 
| 

y=n, 
|| 1 y 9 
2\¢ 6 
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so that P;, ..., S; may be written in the form 


P, = _ oy, +200 | + aC + log~ 
m\ 2 om Ov n 


| 

ly , @ m mfa m 

f + log + — + Bj 
2 2( T log n ( € m i) 


1 fd n 1, .\t 
r n m 


(20) 


nl 


Let the coérdinates of a point with respect to the tetrahedron P,P,P,P, 
of the congruence be x), ¥2, x3, 14, and let the coérdinates of the same point with 
respect to the tetrahedron P,,P.,P,,P,, of the surface S, be 1’, x2", x3’, x4’; 
then 
(21) w(x, Vt + + v,01) + + 139 + , 


where w is a factor of proportionality which may be a function of u and 2. 
Substitution in (21) of the value of y, z, p and o as given in (17) gives an iden- 
tity in 1, 21, p1, and o; and leads to the following equations for the transformation 


of coordinates: 


>, + 
wx, =X — + ) 6 (ave bw X3 


¢ 
WX, = + %3 + OM , 


(29) m 
Ox, = V4. 


The computation for the second focal sheet is quite similar. The differential 


equations which characterize S, are 


“vy “uu — b’z + ( + d 


n nn 


10 
O,=- —G,- | og — 
mLY 2 Ov n° 
¢ lo md 
- 2f + log 
s, 
= 
| 
(23) 
Sun = MNZ+— 
n 
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Make the transformation 


Us = A(U,v) , Vo = p(U,v) , 
where 


in order to make the asymptotic curves of S,, uw. = const. and v2 = const., 


parametric. Then one obtains 


= » = €(Au2us — Purr) 
(24) = — + + — 1 » 


Equations (24) substituted in (23) give the differential equations of S, re- 


ferred to its asymptotic lines. These equations are 


(25) + 2224, + = 0, 
, 
Zum, + + + = 0, 


(26) = Ge 2by = (‘“ 
My € € 


Nu My 
l 
Co = — (b' + 2emn) , co = — — 2emn) , 
de? wy, 
where 
1 fa’ Ny €u l 1 ny 
(27) 
1 
€ 
l 
Ge = Bo 
€ 


The fundamental relative semicovariants of S, are 


Yo = 2, 2 = Zy, + G22, pe = 2, + 


02 = + + + + (b2) us + 


where 
2a2 = + - + 2a2 Fs, 
Au Ay € Ku 
| 
| 
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Using (14) and (28) one obtains 


dod, + boy 
n n n 


V2; 


(29) 
+ + Rope + , 
o=—€ (an, — bop, + + — , 
en 
where 
P; Ay Ku (2) + (b2) us +> 2a2b, 
n Ay Mu 
1A nN, ra) n ly Ny ra) n 
2 Xx E Ou (+ "8 2 + ov (i+ *) 
(30) n 2 \c’ m f’ 
Ay] cn 
= — |  -—G.-- 2 
10 
= 2 Hu 
n 


Let the coérdinates of a point with respect to the tetrahedron P,P,P,P, 
of the congruence be %, x2, 3, %4, and let the coérdinates of this point with 
respect to the tetrahedron P,,P,,P,,P,, of the surface S, be x,, %2, %3, %4; 
then 


where z is a factor of proportionality which may be a function of u and v. If 
the values of y, z, p and o as given in (29) are substituted in (31), there results 
an identity in ye, 22, p2, and o2 which leads to the following relations: 


dod, + b, n 
ti, ( + X2 + P 2x3 a. Nu bo py + 


= “x2 + + 
n 


to 


Ky 
= + Roxs — 


4 2 Ay My x 
Ad 
n 


(38) 
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3. EQUATIONS OF THE EIGHT LINES REFERRED TO THE TETRAHEDRON OF THE 
CONGRUENCE. A STUDY OF THE POSSIBLE COINCIDENCES 


The equations of the directrix of the first kind of P, referred to P,, P., P, 


P,, are* 
(33) 0, x) 2a. 2b, X3 = 0. 
If the values of x,, ..., X4 as given by (22) are substituted in (33), the equations 


of the directrix of the first kind of P, referred to P,P.P,P, are found to be 


xs = 0, 


My 
<0) 


im 2a, 


(34) 


From equations (11) and (12) we find 


+ (tos + (105 9 

(35) 
2( i) 2 ( 


so that the equations of the directrix of the first kind of P,(d) referred to P,P,P,P, 


become 


(36) 0, Lixe Mix3 = 0, 


where 
452 f + log 8m? FG | 
46mLou Gi oy 
4 


The equations of the directrix of the second kind of P, referred to P,,P.,P,,P,, 


| 


(37) 


are 
(38) x2 + On %,=0, = 0), 
2b, <Q, 


and these equations by the use of (11), (22) and (35) give as the equations of the 
directrix of the second kind of P, (d,) referred to P,P.P,P, 


(39) Xo Sing =0, Tiny = 0, 


*E. J. Wilczynski, Projective differential geometry of curved surfaces, these Transac- 
tions, vol. 9 (1908), p. 95, equations (70a) and (700). 


0. 
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where 


. ize) OF .. 
Si = log -§ 3f + log 


- log — 6 — log 


(40) 


In a similar way one may obtain the equations of the directrix of the first kind 
of P, (dz) referred to PyP,P,P,; they are 
(41) vy = 0, Lewy + + Moxy = 0, 


where 


and the equations of the directrix of the second kind of P (d.) referred to P,P .P.P,, 


(42) 


namely 
(433) vi + Sov, = 0, + x4 = 0, 


e'n? _._. 1 oO I. 
Se Bf + log 
| [2.( log ) € Ov 


I, 1 Oo 
| log 4 (/— tog 
len Ou Ges € Ov 


The equations of the axis of P, (4) are* 


he re 


(44) 


(45) v3 = 0, Ay, = 0, 
where 
1 n 
(46) f, + 
d n 


The equations of the ray of ?,(7)) are] 
The equations of the axis of P2(x,) are 


(48) = 0, Kors = 0, 


* i. J. Wilczynski, The general theory of congruences, these Transactions, vol. 16 
(1915), p. 313. 
t Loc. cit., p. 318 


|| 
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where 
and of the ray (re) 
(50) Ve 0), Xs 0. 


Certain relations between the relative invariants Be, Co, Ky, 
7, will be useful in the sequel, and are given here. They may be verified by 
computation : 
Si — mL, = Ky — 2C; 
(51) - nh» Ke 
2B,—-— M, , 
20. — M, . 


The points ?, and ?, will not coincide for all values of w and v, since we are 
excluding the case of congruences with coincident focal sheets.* ‘Therefore the 
tangent plane to S, at P, (x4 Q) is in general distinct from the tangent plane 
to S, at FP. (ts 0). From these two facts one finds that certain pairs of 
lines cannot coincide (this may also be seen from the equations of the various 
lines); in fact, there remain only twelve possible coincidences. Let (dd) 
denote the coincidence of the two lines dy and dy etc.; then the twelve possible 


coincidences and the corresponding analytic conditions are as follows: 


(d,d,): = =Se—M,=0; Le = 7, = Si — Mz = 0; 


(r,%_): Ae = 0; (rex): Ky = 0; 

(d\r,): = M,=0; (d,r2): Le = M2 =0; 
(d,\x.): = M,— Kk, =0; Le =M.—K,=0; 
(d,x,): S; —K,=0; T2 Ss —K,=0; 
(d\r.): =0; (d; 0; 


where the values of Ay, ete., are given in the preceding section. 

The axis and ray of a point on either focal sheet of a congruence correspond 
by duality and the two directrices of such a point also correspond by duality. 
So the coincidences (d,x,), (d,%»), (dir.) and (d,r,) will be called the duals of 
(d,r,), (d,X») and (d,x,), respectively. If a congruence has a coincidence 
then the congruence oblained from the given one by duality has the dual coincidence. 
In fact, the substitution which transforms the coefficients of a congruence into 
the coefficients of the dual congruence transforms the conditions for any coin- 
cidence into the conditions for the dual coincidence. ‘This fact enables one to 


reduce the amount of calculation involved in the study of the possible combina 


* KE. J. Wilezynski, Brussels Paper, pp. 11-12. 


| 
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tions of the twelve coincidences. For example, the canonical form and proper- 
ties of the congruences characterized by the coincidences (d,r,) and (dsr) 
gives by a substitution alone the canonical form and properties of the congruences 
characterized by the coincidences (d; v,) and (d, x»). 

It is not evident that there exists a congruence which has any one of the twelve 
(apparently) possible coincidences, for there might be no solution to the integ- 
rability conditions and the corresponding coincidence conditions taken together. 


However, the congruence which has the coefficients 
(52) m=n=c'’=d=1, a’ =b=c=d'=0, a= — dD’ = const., 


satisfies the conditions for all the coincidences, so that there surely exists at 
least this one solution for any group of equations obtained by using any com- 
bination of the twelve coincidences. ‘The congruences characterized by (52) 
will be referred to as coincidence congruences in the remainder of this article. 
Since m ~ 0, n ~ 0, neither of the focal sheets of a coincidence congruence de- 
generates into a curve, and since c’ # 0,d ¥ 0 neither of its focal sheets is de- 
velopable.* We shall throughout consider only such (mnc’d #0) non-degenerate 
congruences. ‘The coincidence congruences have four (d;, di, dy, d») of their 
eight lines indeterminate, so, that the coincidences in this case have but little 
meaning. We shall show that certain sets of coincidences give rise to congruences 
for which the eight lines are determinate, while for certain other sets some of 


the eight lines are of necessity indeterminate. 


!. ‘THE PROPERTIES AND CANONICAL FORM OF CONGRUENCES POSSESSING COIN- 
CIDENCES (d)7;) AND 


The conditions which give these coincidences are L; = M,; = S; = 7; = 0. 


from these equations and (51), we find 


53) B,=0, K,=2C,, Fy =G, = 
Using (53) in L; = 0 and M, = 0 gives 


ra) 
ov\ 1 
whence by integration 


(54) C; = V(v), Che = U(u), 
m 


* E. J. Wilczynski, Brussels Paper, p. 28. | 


= 0), — 
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where \(v) is an arbitrary function of v alone, and (’(#) is‘an arbitrary function 
of « alone.* 

Integration of C; = \V(v) and 4, = 0, using (11), gives 


(55) = V;(v)U,\(u), 


where V’;(v) depends upon \(v). From (54) and (55), one finds that 


UV." 
so that /,, 0. Consider the most general transformation (3) which leaves 
the set (1) invariant; since (c’ ) , we can choose dA, yw, @a,, 8, and A 
in such a way as to make 
(00) e=], 


and this relation will be left unaltered by the subgroup of (3) for which \ = 


C2 
—, = —, where ¢ and are arbitrary constants. 
Vay B 
From (55) and (56) we find that om y —— ; since =- , we can 
m V2 m ay, 
make d = m, and thereafter use the subgroup of (3) for which a, = 6, = c, 


AX = co, w@ = C3 where ¢, cz and c; are arbitrary constants. 


We now have, using the fifth integrability condition, 
(57) f=0, d=m, =n. 


These values substituted in (54) give mC, = g(u); moreover we find by calcula- 


tion that = —, and from = Othat m = so that g(u) = 
2m 


m, while 


wo 


l, 


2(lv + hh) 


m = yY(v). Thenm = lv + 1, and Ci 


*In the remainder of this article, arbitrary functions of u alone or of v alone are indicated 
by this same notation. 

t Cf. E. J. Wilczynski, Brussels Paper, pp. 42-3. Similar transformations will often be 
made in the remainder of this article; in every case, the sub-group of (3) which is available for 


further transformations will be indicated. 


| 
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There are two cases to consider. 
First, 12 = 0; then m = const. and we find that the transformations at our 
disposal enable us to make 


(58) d=m=c'=n= 1. 
Using (58) the integrability conditions give 
(59) a’=b=0, c=d’'=0, b= —a= const., 


and (57), (58) and (59) characterize the coincidence congruences. 
Second, lz # 0; now (m) = m, so we can make /s = 1, and thereafter 
C\C2 


use the subgroup of (3) for which 


Nowd =m=v24+1, 2C, = =K,= log d whence c’ =n = U3(u). The 
v+ il; Ov n 


sixth and seventh integrability conditions now become a + b’ = O and n 


uu 
l 

+ n,, = 0 sothatn = c’ = but since (c’) = — c’ we can make]; = 


Using 
(60) f=c=d'=0, d=m=vth, c’=n=u+h 


in the integrability conditions, we obtain 


(61) a’=b=-—1, b’=-a, 
where a, = 2(v + 1), a, = —2(u + 12), so that 
(62) a=v"— 4. 2(Lyw Isat) + ls 


where /;, /, and /; are arbitrary constants. 
If we use the values given by (60), (61) and (62), we obtain lz; = Mz = Sy 


= 7, = 0, Ki = ——-, K, = ——\~,,so that these congruences possess the coin- 
u + 


*Ifl; = 0, this cannot be done; but as in the first case above, we can show that if /; = 0 
the congruence is a coincidence congruence. 


a, = B, A = C2, = CiCe. 
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cidences (d,d,), (d.4,), (d,r.) and (d,r,) in addition to the given pair. To 

recapitulate: Jf a non-degenerate congruence possesses the coincidences (d,r,) and 


else it is a congruence whose canonical form is 


(63) f=c=d'=0, d=m=vth, cl=n=uth, 
a=b=-—-1, =v? — 4 — hu) + i; 


. . . 
where l;, lz and ls are arbitrary constants, and possesses also the coincidences (d,d.), 
” ” 
(d.d,), (d,r.) and (d,r,). 


‘The congruences dual to (63) have the canonical form* 


(64) f=c=d'=0, d=m=ut+h, 
b = a’ = 0, b’ = y* — 2(lw Ist) + ls, 


and possess the coincidences (d3x,), (dix,) and (d5x,). 

Since there exist congruences, other than the coincidence congruences, which 
possess either one of the sets of coincidences noted above, there exist congruences 
which possess any one of the coincidences separately. We will consider in the 
next section the congruences which possess the coincidences (7:72) and (72%) 
and we shall show that there exist congruences, other than the coincidence 
congruences, which possess that pair of coincidences. Using that fact, we state 
that there exist congruences, besides the coincidence congruences, which possess 
any one of the twelve possible coincidences. 

The coefficients given in (63) make kK; = sedi ~ 0; if we seek those con- 

1 

gruences which possess the coincidence (72x%,), which makes K, = 0, in addition 
to the coincidences (djr,) and (d‘7,), we find that they are coincidence congru- 
ences. More generally the coincidence congruences are the only congruences 
which possess all of the twelve coincidences. ‘This justifies the term “coincidence 
congruences’’ which we have chosen for them. 


5. ‘THE PROPERTIES AND CANONICAL FORM OF A CONGRUENCE POSSESSING COIN- 
CIDENCES (71%2) AND (2172) 


Let us consider now a congruence which has the coincidences (7:72) and (x;72). 
By integration of the associated conditions 


=(0 =( 
, ), 


n 


Su 


c 


* E. J. Wilczynski, Brussels Paper, pp. 24-28. 


(d,7.) it ts either a coincidence congruence which possesses all twelve coincidences, or 
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one obtains 


Transformations (3) enable us by a proper choice of A, uw, @,, 8, to make 


V(v) = 1, U(w) = 1. We may use thereafter the sub-group of (3) given by 
l 
VQ, 


We now have 


(65) c’=me’, d=ne 


Here, as previously, we exclude congruences for which muc’d = 0. 


From (65) we find 


f ‘ (2) 
Cun = — 2m, fy — + 
d., = (n, — nf.)e~”, = (n., — 2n,f, — + 


If these values are substituted in the sixth and seventh integrability conditions, 


we obtain the equations 


— myfy — ma = — e — nf, — nb’), 


— ma = — e’(n,, — n,f, — nb’); 


un 


un 
whence follows either f = 0, which means that the congruence is all’-congruence, * 
or else 


Myy — Myf, —-ma=0, 


(66) 
Nw — Nf, —nb’ =0 


’ 


we shall call such congruences harmonic. 7 
If f = 0, equations (65) reduce to 


(67) c'=m, d=n. 


*E. J. Wilczynski, Brussels Paper, p. 46. 
+ The nets of curves u = const. and v = const. on Sy and S; are in this case harmonic 


conjugate nets. See E. J. Wilczynski, Geometrical significance of isothermal conjugacy of a net 
of curves, AmMericanJournalof Mathematics, vol. 42 (1920), p. 215. 


2 —f d f 
=Vivje? , —=U(u)e 
m n 
] 
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From (67) and the integrability conditions it follows that 


403 
f=c=d'=0, b=—un,, a’= My, 

(GS) Mun Noy = ma+ nb’, 
a, = 3m,n + mn,, 


3mn 


The canonical form given here is the same as the canonical form for identically 
self dual congruences, or for congruences which belong to linear complexes. 
Using (66), we find 


c=f,, d’=f, b=—ne', a’ =—mge’, mn(l—e”) =f, 
n n nn n 


The results of this section may be stated as follows: 


If a non-degenerate congruence possesses the coincidences (1x2) and (xr2), 


ii 1s 
etther a W-congruence whose canonical form is given by (67) and (6S) in which case 
it ts identically self-dual and belongs to a linear compiex; or else it is an harmoni 


congruence whose canonical form ts given by (65), (66) and (69). 
0. 


THE PROPERTIES AND CANONICAL FORM OF CONGRUENCES WHICH POSSESS 
COINCIDENCES (d,d,) AND (d,d.) 


We consider finally the congruences which have the two coincidences (did5) and 
(dd.). ‘The conditions to be satisfied are 
(70) M,=90, L,=T,=35,—M.=0. 
Using (51), equations (70) may be replaced by 
(71) 0 + log 6°m? F\G 0 
+85 og 6°m? F\G, } = 0 ; 
Ou 5 Gy ov f 
(72) Lo => 


log 


(73) M,=2h,, f + log bm? +6 


0 of, log €°n? 0: 


* E. J. Wilczynski, Brussels Paper, p. 28 and p. 43. 


| 
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74) += tog 2 (25 — tog ; 
(74) Me= 2C2, f + log — en 


(75) Si = 5 On log G, 3f + log = = 2f — log ; 


(76 oR en? 1o F, of lor 3%! 
If we multiply (75) by 6 and subtract (71) from the result, we find 
log 6°e? = 0 
og = 0. 


In a similar fashion from (72) and (76) one obtains 


log = 0 , 
so that 
(77) be = k, where k is an arbitrary constant. 


Using this relation, the six conditions reduce to the four given by (71), (72) and 


the two following: 


‘ ou \- og om de og 

(79) log f — log en*F.G2}) =0. 
Adding (71) and (78) gives, using (11), 

while subtracting (7S) from (71) gives 
(S1) F iy — 26F\G, + + 6F,, = 0. 


Adding (SO) and (S1) gives 


+ (Fi — Gi) 6, +8 > = (), 


| 
| 
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») 
which, since + G; = —2C), =- B,, becomes 
6 


Cin + Bi, = 0. 
But from (11), 


so that 
2) f+log6}) =0 
Oudv\’ 


By a similar calculation, from (72) and (79), one obtains 


83 + loge} = 0. 
Oudv f 
From (S82), (S3) and (77) it follows that /,,, = 0, whence 
=1, andk= +1. 
mn 
But « = += 1/6 enters in such a way into equations (71), (72), (78) and (79) 
that the same set of equations results from using k = —1 as from using k = 1. 
So there is no loss of generality in assuming k = 1. Use the transformation 


group (3) to make 
(84) f=const., c=d’=0; 


these relations will be left unaltered by the sub-group of (3) given by 


Ci Ce 
Vay VB, 


where ¢; and cz are arbitrary constants. 
From (84) and (S82), it follows that 


P= — — = — ; 


then by using transformations (3) we can choose @ and £8 so as to make 


(85) d=m, 


— 
| By = ——(f-+ loge 

( lu + it Oudv ft+ og 0 

| 
| 
| 
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and thereafter use the sub-group of (3) given by 


where and are arbitrary cons‘ants. 
Since mn — c'd = f,, = 0, we find further 


(86) c= 


These equations show that the developables of the congruence intercept isother- 


mally conjugate nets on both sheets of the focal surface. * 
Subtracting (80) from (81), and using (11) and (S5), one finds 


(87) Mun + My = 0. 
By a similar computation based on the two equations corresponding to (80) 
and (81) one obtains 


(SS) Nyy + Ny = 0. 


The integrability conditions become 


(89) b=—m,, a’=-—n,, b’=-—a, 
(90 a, mn, A, = — MN. 

u ov 


If we differentiate the second equation of (90) with respect to v, and the first 
equation with respect to u, and subtract one resulting equation from the other, 
we obtain 

(91) (mn) un + (Mn), = 0. 


If we perform the differentiations indicated in (91), and use (87) and (88), we 
obtain 
My mM, 


(92) =— = p, 
n Ny 


* E. J. Wilczynski, General theory of congruences, these Transactions, vol. 16 (1915), 


C1 Ce 2 
% a, = B, = C3 
3 3 
| 
p. 322 
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where p is a function of “ and v which we may assume to be real, if the focal 
sheets of the congruence are real. From (92) we find 


Muy = + PU yy, Moy = — — » 
whence, by (87), 
(93) PuNy — = 0. 
Again we find 
Muy = Py Ny + pny, = — — 


whence, by (88), 
(94) + = 


Equations (93) and (94) are two homogeneous equations which , and 1, 
must satisfy. The determinant of the coefficients must be zero;* i.e., p, + 


p, = 0, whence p, = p, = Osince pis real, and p = k, wherek is an arbitrary 
real constant. Then k = —m,/n, = —b/a’, using (2); so by using (3) we can 


make k = 1 (k = 0 gives the coincidence congruences) and thereafter use the 
Subgroup of (3) for which 


Equation (92) becomes 

(95) mM, = —Ny, M, = N,; 

using these values, equations (90) can be integrated and give 
(96) a= n? 


where & is an arbitrary constant. 

The congruences which possess the coincidences (d,d,) and (d/d;) are 
identical with the congruences which arise from the Riemann sphere representa- 
tion of functions of a complex vairable according to one of the methods discussed 


* Or else nu = n, = 0; then = const. and if m # 0 we can make m = 1 and are led to the 
coincidence congruences, which form a special case of the general set of congruences possessing 


the given coincidences. 


“ore 

| 
| 
| 
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by Wilczynski.* We shall call them bidirectrix congruences since the four 

directrices associated in general with a line of the congruence in this case 

reduce to two. Thus the property of possessing coincidences (d,d;) and (d\d;) 

characterizes completely this set of bidirectrix congruences, discovered by Wilczynski. 

Using Wilczynski’s results in connection with the theorem just stated we see 

that a non-degenerate congruence in the real and distinct focal sheets, which possesses 

the coincidences (d\d,) and (d\d;), has also the following properties: 

(a) it ts a W-congruence; 

(b) tts developables intercept isothermally conjugate nets on both sheets of the focal 
surface; 

(c) the asymptotic curves on both focal sheets belong to linear complexes; 

(d) corresponding asymptotic curves belong to the same complex ;t 

(e) the directrix curves of the two sheets of the focal surface correspond to each 
other: 

(f) on each sheet of the focal surface, each of its axis curve tangents 1s conjugate to 


one of the ray curve tangents. 


Moreover 7t has the canonical form 
c=d’'=0, d=m, c' =n, a’ = —b=—n,, a= — db’ = n? — 


where 


and k is an arbitrary constant. 


* Line geometric representations for functions of a complex variable, these Transac- 
tions, vol. 20 (1919), pp. 283-298. 

A set of properties characteristic of a class of congruences connected with the theory of func- 
tions, these Transactions, vol. 21 (1920), pp. 409-445. 

+ This property is added by the author and may be proved as follows. Denote byT the 
linear complex to which the asymptotic curve v; = const. on Sy, belongs so that the linear com- 
plex osculating 1; = const. at any point is’. In view of the coincidence property, the two 
linear complexes osculating u, = const. and v7; = const. which pass through P, have in common 
a linear congruence whose directrices are the directrices of the first kind of Py and P:. All the 
lines which meet these two directrices belong to , so that [ may be described as follows: 
consider the directrices of the first kind of the points on an asymptotic curve on one focal sheet, 
and the corresponding directrices of the first kind along the corresponding asymptotic curve on 
the other focal sheet; then the linear complex to which the first asymptotic curve belongs is 
made up of the totality of lines which meet a pair of corresponding directrices. Evidently 
the corresponding asymptotic curve on the other sheet also belongs to the same linear complex. 
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ASYMPTOTIC PLANETOIDS* 


BY 
DANIEL BUCHANAN 


1. Introduction. Near the vertices of the equilateral triangles described in 
the plane of Jupiter’s orbit and on the line joining the Sun and Jupiter as base 
are to be found six of the planetoids. Four of these,t Achilles (588), Patrocles 
(617), Hector (624), and Nestor (689), oscillate about one vertex in much the 
same way as the mythological luminaries of the same names circulated about the 
walls of Troy. Two planetoids are found at the other vertex, but their names 
are unknown to the author.{ These six planetoids are in the vicinity of two of 
the five well known points of libration in the problem of three bodies.§ The 
three other points lie on the line joining the Sun and Jupiter, one point lying be- 
tween the Sun and Jupiter, another on the side of the sun remote from Jupiter, 
and the third on the side of Jupiter remote from the Sun. One of these straight- 
line points in the case of the Sun and Earth, viz., the one on the side of the Earth 
remote from the Sun, has a physical significance in that it may account for the 
“Gegenshein.”|| ‘The equilateral triangle points of libration were considered by 
Lagrange in his celebrated prize memoir of 1772 as ‘‘pure curiosities,’ but re- 
cent astronomical discoveries show that these points likewise have some phys- 
ical significance attached to them. 

The object of this article is to determine orbits for these planetoids of the Sun 
and Jupiter which will approach the equilateral triangle points of libration as the 
time approaches infinity. As these orbits are asymptotic to the above men- 
tioned points of libration, we have designated the planetoids which move in such 
orbits as ‘‘asymptotic planetoids.”” These planetoids are considered to be of 
appreciable mass, and their perturbations upon the Sun and Jupiter are deter- 
mined. ‘They are assumed to move in the plane of Jupiter’s orbit. 


* Presented to the Society, December 31, 1919. 

Tt Marcolongo, J/ problema dei tre corpi, Scientia, vol. 1 (1919), No. 8. 

ft I am indebted to Professor E. W. Brown for this information, as he communicated it 
to me when this paper was presented to the Society. 

§ Lagrange, Collected Works, vol. 6, pp. 229-324; Tisserand, Mécanique Céleste, vol. 1, 
chap. 8; Moulton, Celestial Mechanics, pp. 309-318. 

|| Gyldén, Sur un Cas Particulier du Probleme Astronomique, vol. 1; Moulton, Celestial 


Mechanics, p. 305. 
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In plain mathematical language, divorced from every astronomical applica- 
tion, the paper treats of two-dimensional asymptotic orbits near the equilateral 
triangle equilibrium points in the problem of three finite bodies. The paper con- 
cludes with numerical examples of orbits in which the ratios of the masses are 
not those of the Sun, Jupiter and the planetoids, but are chosen so that the or- 
bits near each point can be drawn to the same scale. 

Several classes of asymptotic orbits have already been obtained, but with one 
exception, number (6) below, they belong to the particular case of the problem 
of three bodies in which one body is infinitesimal and the finite bodies move in 
circles. The following is a list of the solutions which have been ob- 
tained 

(1) Warren determined two-dimensional orbits which are asymptotic to the 
points of libration lying on the straight line joining the two finite bodies. 
(American Journal of Mathematics, vol. 38, No. 3, pp. 221-248.) 

The remaining cases were found by the author of the present paper. 

(2) Two- and three-dimensional orbits which are respectively asymptotic to 
the two- and three-dimensional periodic orbits near the straight line equilibrium 
points determined by Moulton* in the chapter on Oscillating satellites. (A me r- 
ican Journal of Mathematics, vol. 41, No. 2, pp. 79-110.) 

(3) Two-dimensional orbits which approach the equilateral triangle points of 
libration. (Transactions of the Cambridge Philosophical 
Society, vol. 22, No. 15, pp. 309-340.) 

(4) Three-dimensional orbits which are asymptotic to the three-dimensional 
periodic oscillations near the equilateral triangle equilibrium points which 
were determined by Buck.j ‘This paper forms part of the paper mentioned 
in (3). 

(5) Three-dimensional orbits which are asymptotic to the isosceles triangle 
solutions{ of the problem of three bodies. (Proceedings of the Lon- 
don Mathematical Society, ser. 2, vol. 17, No. 1, pp. 54-74.) 

(6) Two-dimensional orbits which are asymptotic to the straight line equi- 
librium points when the three masses are finite.§ 

2. The differential equations. Let m1, m2, and m; represent the masses of the 
three bodies and let M denote their sum. Let a system of rectangular axes be 
chosen having the origin at the center of mass of the three bodies and the plane 
of their motion as the plane of reference. Let the axes rotate about the origin 
with the uniform angular velocity and let the coérdinates of the three bodies 
m, be denoted by (x;, y;),7 = 1, 2,3. If the bodies are subject to their mutual 


* Periodic Orbits, chap. 5. 

+t Buck, chap. 9 of Moulton’s Periodic Orbits. 

ft Buchanan, chap. 10 of Moulton’s Periodic Orbits. 
§$Rendiconti del Circolo Matematico di Palermo, vol. 45, pp. 1-25 
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attractions according to the Newtonian law, then the differential equations 


which define their motion are 


” 9 9 1 
— — n*x; = —, 
m; OX, 
1 aU ‘ 
vi + — ny, = — — = I, 2, 3), 
(1) m; OY; 
>» [mime moms m3m, 
Ti2 723 31 
ry = — 4? + Oi — %)?) = 1,2,3; #4), 


where k? is the factor of proportionality and the accents denote derivation 
with respect to f. 
These equations admit the vis viva integral 


3 

(2) + — x37 — v3) = 2U+-C, 

where C is the constant of integration. 

It is shown in treatises on celestial mechanics* that the equilateral triangle 
configuration with proper initial components of velocity is a particular solution 
of the differential equations of motion (1). If the units of distance and time are 
so chosen that the mutual distances and k? are unity, then the angular velocity of 


rotation 7 must satisfy the condition 
= M. 
Since the unit of mass is so far arbitrary, it is possible, without loss of generality, 
to put M = 1. Hence n? is also unity. 
If the axes are initially orientated so that the x-axis is parallel to the side of 
the equilateral triangle joining m2 and m3, then the values of the coérdinates for 


the particular solutions are either 
Configuration I (Fig. I), 


1 = (m3 — me), = (m2 + ms), 


l l , 
(3) = 5m + mz, w= 


=- (5 my + me }, 


* See Moulton’s Jntroduction to Celestial Mechanics, pp. 309-11 
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or 
Configuration IT (Fig. 2), 
1 /; 
= (m3 — me), m= V3 (me + M3), 
] | 
V3 m. 
x 
X / 
m / 
Fic. 1.—Configuration I. Fis. 2.—Configuration IT. 


Since the origin of coérdinates is at the center of mass, it follows that 
+ + M3X3 = V, + Moye + M3V3 = O, 


and therefore the coérdinates of one of the bodies may be eliminated from the 


differential equations. Let us suppose x; and y; are thus eliminated by 


] 
(5) + = — — (myn + Moye). 
Wa Nz 


When these substitutions are made in (1), and k? and » are both put equal to 


unity, the differential equations of motion become 
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M(x, — Xe) (992, + + 
Pie Piz 
” M2(V; — V2) (331) 
P| 
— X1) + (ine + 
Pi2 Po 
(6) ” — V1) + (Me + M3) V2 
Pie Pox 
” 1 2 
pio = — Xe)? + (11 — 
l oil /2 
pis = — + + — + 
M13 


= + (1 — + + my) 
We 


3. The equations of variations. Let 


v1 = 5 (ms — mz) + € i, V3 (ite + m3) + €%, 
(7) 
= my 4. + le, 2 5 Vom +. € V2, 
where 7,. . ., v2 are new dependent variables, and ¢ is an arbitrary parameter. 


The additive quantities e7),. . ., ev, denote the components of displacement of m1, 
and mm». from the vertices of the equilateral triangle in configuration I (Fig. 1). 
We shall not consider configuration II (Fig. 2) in detail as the coérdinates in 
Fig. 2 differ from the corresponding ones in Fig. 1 only in the sign of x/3. Since 
4/3 does not occur explicitly in the differential equations, we may change the sign 
of V3 in any solutions which we obtain for configuration I and thereby obtain 
the corresponding solutions for configuration IT. 

When equations (7) are substituted in (6) the right members of the differen- 
tial equations can be expanded as power series in e which will converge for e 
sufficiently small. Since the expressions in (7) that are independent of ¢€ are 
particular solutions of equations (6), there will be no terms independent of ¢ after 
the substitutions (7) have been made. Hence the factor «can be divided out of 
the equations, and when m; is replaced by 1 — m -- my the differential equa- 


tions (6) become 


| 
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Dz — uy — (1 — + Ou2 — mare 
4 4 f 2 

= UP + +... +f + ..., 


9 


= Vf3 (eV + +... + .--], 


9 3+/3 9 30/3 
— — my, + G —3+- m) Us — (2p m:) Ve 
3V3 9 3V3 
mu, + + (20 m) me + +m) Ve 
= + +... + + ..-], 


(S) 


where D denotes the operator d/dt, and ut... .opVeq (& = 2,...,m,...) are 
homogeneous polynomials in 7, us, %, v of degree k. In all these poly- 
nomials, +/3 occurs only as a factor of the odd powers of v; and z, considered 
together. 

If we consider only the linear terms in (8) we obtain the equations of variation. 


They are 


uy — 2D + ~ = (1 — + One — ve = 0, 
2D — —— (1 — 2me) 7 u D? — —)}) % — —— motte + Ove = 0, 
(9) 
74 — — 3 +5 — (2p m:) te = 0, 


-~ + j my + (2 + ~ ; m) Uz + = 


4. The solutions of the equations of variations. ‘The solutions of the equa- 
tions of variation can usually be obtained by differentiating the generating solu- 
tions, equations (3), with respect to the arbitrary constants which denote the 
initial time and the scale factor.* This is not possible in the case under con- 
sideration since the generating solutions are constants. 

Equations (9) are linear differential equations with constant coefficients, 
and in order to obtain their solutions we consider the operator D as an alge- 


* Poincaré, Les Méthodes Nouvelles de la Mécanique Céleste, vol. 1, chap. 4. 


1922] ASYMPTOTIC PLANETOIDS $15 
braic quantity and equate to zero the determinant formed from the coefficients 


of 1, Us, 21, v2 in (9). This gives 


— — 2) (1 — 22), 0, 
(10) 
9 
2D ¥ (1 — 2m), Dt — ~ 0 
A= 
= () 
9 34/3 34/3 
34/3 3/3 , 


This determinant reduces to 


A = D? | + + 1 + 


)7 


af 
+ + + mam) = (). 


+ 4+ D? 
(11) 


If we neglect the factor )*, the preceding equation is a cubic in /)* and its dis- 
criminant is 


+ + mg)? [27 + + m3) — 1) 
) 


Since 1, + ut2 + m3 = 1 and since m1, ms, and ™s3 are all positive, this dis- 
criminant is positive, and therefore the cubic has one negative root and two con- 
jugate complex roots. Let these roots be Xj, A;, and A; where 


Then the solutions of (11) are 
(13) D = 0, + de, + Az, 


and the solutions of the differential equations (9) are therefore 


Aut dot Nat — 
Uy = ay + aot + + aye + ase” + ave + + ase 
(0) (0) Al (1) —Ail (2) (2) 
= a1 + By act + By’ ase™ — + Base” — By ace 
(3) (3) Ast 
+ B, ac” — a 
(0) (0) (1 Aut (1) —Ail (2) Rot 2) - 
(14) we = a1 + Yo aot + 7; age + Yo ae” + ase” + ae 
(3) Aal (3) —Asl 
+ y; + Yo ase 
(0) s(1) —Ail (2) Aol ¢(2) — Act 
= 5," aot 0; ase" 05 ase 0, — 05 age 
s(3) Ast 
+ 6; ae" — 65° ase 
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where a, ..., as are the constants of integration. ‘The remaining symbols 8,’, 
...,6, denote constants which are defined by the following equations and prop- 


erties: 
0 l Ss ‘3 2 
B; (l — m,) — (mm, — V5 me) 
2 
+3) + 4inems} |, 
Yim, — 1) 
(at, — 1) (V3 + + V/3 (me + m3) + 
— 
l S(m3 — moe) /3 
- (1 — + : — + 
2 Q(m, — 1) 3 
— Moe my my, + 
V — 1) — V 3(m — 1) 
3) 373 [(m3 — me) + 2y 321 — m,)A, + — ms) | 
— — Ym) + 6V/3(m3 — me)dA; + 27(1 — m) 


9 3 — 2m, 


Qin Qin» 
> 3 3 1 — 2m F 

) A; + 2 By (¢ = I, 2, 3). 


‘The following pairs of constants differ only in the sign of 7/3: 


‘The following pairs are conjugate complex: 
I 2 (3). 2 3). 3), ‘ 


The solutions for ) in equations (13) are called the characteristic exponents. 
Since the original differential equations (1) admit the integral (2), then two of 
the characteristic exponents will be zero.* It was to be expected, therefore, 
that two roots of (11) should be zero. Poincaré+ has also shown that if equa- 


* Poincaré, loc. cit., p. 188. 
t Loc. cit., p. 69 


B,’, Bs’; 3; (¢ = 1,2. 3). 
| 
| 
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tions (1) do not contain ¢ explicitly, which is always the case in all problems of 
mechanics in which there is a conservative system, then the characteristic expo- 
nents are always equal in pairs but opposite in sign. It is this property of the 
solutions that makes the determinant (11) even in ). 

5. Construction of Asymptotic Solutions in ¢ ““ and ce“. We shall now 
construct solutions of the differential equations (4) which are asymptotic in the 


sense of Poincaré,* that is, each term of the solution must have the form 
P(i), 


where c is a constant having its real part different from zero, and P(t) is a periodic 
function of ¢ or, in particular, a constant. Such solutions will therefore approach 
zero as t approaches + “ or — % according as the real part of c is negative or 
positive, respectively. 

The only terms of the solutions of the equations of variation which are 
asymptotic in the sense just defined are those ine~“ ande~’”. ‘These terms 
approach zero as t approaches — “ or + « according as the + or — signs are 
taken with and A3. 

In this section we shall construct the solutions which approach zero as t ap- 
proaches + %. ‘These solutions will obviously involve powers of e “” and 
e ’. Inthe next section we shall show how the solutions in e~™ and e**™ 


. — al — Asl 
can be obtained directly from the solutions in e and e 


Only the formal construction of the solutions is considered in $§ 5 and 6, but 
the convergence of the solutions obtained is established in § 7. 

We propose to integrate equations (8) as power series in ¢ and this is the reason 
why ¢ was introduced in equations (7). Accordingly we substitute 


L L 
(15) = > v;; €? = 1, 2), 


in equations ($8), and as we shall have frequent occasion to refer to the resulting 
equations we shall cite them as (S8’). These equations (S’) are to be satisfied 
identically in « and we may therefore equate the coefficients of the same powers 
of «. In this way we obtain sequences of differential equations in ;; and 7,; 
which can be integrated step by step, as we shall show, subject to suitable initial 
conditions and restrictions. 

The first two steps of the integration will be considered in detail, and then an 
induction to the general term will be made to show that the process of integra- 
tion can be carried on indefinitely. 


* Loc. cit., p. 340. 


7=0 j=0 
| 
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Step 1. Terms in (8) independent of The termsin (8’) which are indepen- 
dent of ¢ are obviously the same as the equations of variation (9) if the subscripts 
on the dependent variables in (9) are altered so as to read uy, V19, Meo and Veo 
instead of %, “2, 11, %2. Hence the solutions of these equations are the same as 
(14) with the corresponding changes of subscripts. Now only two of the ex- 

— at —Asl 
ponentials, viz., ¢ and e “, approach zero as t approaches + % and we 
therefore put equal to zero the arbitrary constants associated with the other 


exponentials. The solutions at this step are, then, 


(0) (0) 2) — (; (0) 
) 
(2) _ (0) Aol (3) _(0 —A3l (2 0) dot 3 (0 — ot 
= Yo e + Y2 Gs - 6. e — 6,” ag 
(0) (0) . (2) ~ 
where a, anda, arearbitrary, and §,,...,6, are defined as in equations (14). 


It is evident that the complementary functions at all the succeeding steps of 
the integrations will be the same as (14) and that after the exponentials are re- 
jected which do not have the proper form, two constants of integration will re- 
main undetermined. Hence, besides rejecting all the exponentials except e~™ 
and e ™, it is necessary to impose two initial conditions upon these solutions 
so as to determine the arbitrary constants arising at each step of the integration. 


Let us suppose that 
(17) mu, (O) = (0) = 
where «, and y are arbitrary. Since #; and #2 are multiplied by the arbitrary 
parameter € in (7) we may put « or y equal to unity. Let us suppose « = 1. 


When the conditions (17) are imposed on (15) it follows that 


(0) = 1, u;; (0) = 0, 
(18) 1 lj ) 


| 
S 

= 


(0) = uo; (0) 


Now applying these conditions to the solutions (16) we obtain 


2) 
a - a, = 
, (3) (2) 3) 


) 


Since y,” and y;” are conjugate imaginaries, it follows that a,” and a,” are like- 
wise conjugates. 


In order to unify the notation, we put 


(9) 


(0) (10) 2) 10) 2 (2) (10) 
a, = Aig — Bo ae = By, Yo % = Ay, — a = By , 


(19) 


(9) (o1) (3) .(0) — (01) ( s(3) (0) _ 01) 
Bs Gs = Bin > Yo =n, = Bay 
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3) (2) (3 2) 


Since a,”, ay; By’, Bo’; Yo» Yo 3 52» 5, are conjugate pairs, it follows 
that the symbols in (19) which differ only in a permutation of superscripts are 
likewise conjugate pairs. This notation will be adopted in the sequel to 
denote conjugate complexes. Thus at and —" (j # k), are conjugate 
imaginaries. 

On employing the above notation, we obtain for the desired solutions at this 


step 

(20) u;, = + » Vio = B; + B; e ig = 2), 
10) (01) 

where a;,_, .,B;, are linear in ‘y. 


If we put 
w= 


and suppose that 
= fai? + 1 Bi? = + 1 


then the solutions (20) become 


— pl 1) 1) ° 
= [aj, cos vt + bj,’ sin vt], 


(21) 
=e cos vt + sin vt]. 


This second form of the solutions is more convenient than (20) for numerical 
computation but is more cumbersome in obtaining the solutions at the succeed- 
ing steps. 

Step 2. Coefficients of € to the first degree. When the solutions (20) have been 
substituted in (S8’), the differential equations obtained by equating the coeffi- 


cients of € to the first degree are 


3 3 30/3 

(v: ) un — 2D + (1 — + 0 uy — = Un, 
3/3 9 34/3 


(2 m) = Us, 


37/3 9 30/3 
+ + (2 + Us, + (0: m) Vo = Vo; 


= 
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20) — 2r2z —(A2 + (02 
Un = 4 ¢ 92st 
9 af = {he (02) .—2Asi 
Vin Bee e Aol + e t A3) Bi e 


(23) 
(20) +A 02 —2dAz 
l 1 > As e + As.” e we + AS 
, (20) —(Az + As) (92) 
where A,;", ..., By, are quadratic expressions in y having constant coefficients 


The constants like A,; and A,,” which have the same subscripts but which have 
their superscripts reversed are conjugate imaginaries. The other constants 

11) 

The complementary functions of equations (22) are the same as (14), and the 


»(11) 
..+, Bg, are real. 


particular integrals can be found by the method of the variation of parameters 


as at the previous step. These particular integrals are 


= + kd3)t (7 + k 2, ; l, 2), 
ji,k=0 
(24) 
j, 0 
where 
k) (jk) 
As 
(as) Al[— + ] A[— + 
20 


The preceding symbols are defined as follows: 

A[ — (jde + kds)] denotes the determinant A in (10) when D has been re- 
placed by — + FAs); 

i = 1, 2, 3, 4, denotes the preceding determinant A[ — (jd. + kds;) ] 


when the elements of the 7th row, reading from top to bottom, have been replaced 


jk) (jk) jk) (jk ° 
by Ay, Bn» An » Ba » respectively. 
lhe constants in (25), viz., a, , ..., 62; are found to have the same form 


as the constants in (23). Thus —" and as) are conjugate complex if 7 ¥ k, 
but real if 7 = k. 

If the right members in (22) contained any constant terms, or terms in ¢ or in 
=A! F = 1, 2,3, that is, terms which are the same functions of ¢ as the terms in 
the complementary functions, then the denominators in (25) would vanish and 
the particular integrals would not take the form (24), but would contain terms 
in ¢, 2? ort e***. Hence, such terms would appear in the particular integrals 
if and only if the right members contained constants, terms in ¢ or ei re- 
spectively. In order, then, to show that there are no terms in the particular 


integrals which do not satisfy Poincaré’s definition of asymptotic solutions, it 
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is sufficient to show that the right members do not contain any terms which 
are the same functions of ¢ as any part of the complementary functions. 

To obtain the complete solutions of (22), we combine the complementary 
functions, equations similar to (14), and the particular integrals (24). The par- 


ticular integrals have the desired form for asymptotic solutions, but we must re- 
ject all the terms of the complementary functions except those in e “ and 

~*_ "This is possible, of course, by equating to zero all the arbitrary constants 
of integration except those associated with the exponentials e~ ande~™”. The 


desired solutions are therefore 


+ 
Ho, = 4. ak” e »— Ast +a ¢ (a2 4+ 
= — at” = 5,” at” Asal 20) e Bo» + 
+ 


where a, and a, are the constants of integration. By virtue of the initial 
conditions (18), these constants must satisfy the equations 


1) (1) (20) (1) 
ag ta, = — lay +a + = 
‘ 2 ‘ 
a) vs” al! — fay? 20) +a ast?) 
or 
(1) ys 1) (1) (1) (2) 
ay — RY? — RY ay, RY — RM x 
2 


(20 (02) (20) (02) ) 


) 2 2 2 . . (11 
Since a,; , @, and , are conjugate pairs and a,, 


(1) 1) . . (2) (3) 
R\” and are likewise real. But as and 
(1) (1) . 
that a, and a,’ are also conjugates. 


In order to unify the notation, as at step 1, we put 


(11) 
and a,, are real, 
are conjugates, it follows 


(2) — fio) (2) — 2) — (10) 
(1) — (01) (3) — (01) (2) 4) — (01) 
a, = > — ag = Bi’; Yo = — 6° = Boy 


and the desired solutions at this step take the form 


9 


Gk) kys)t 
“;, = af” + 


(26) re 


jk) ,— (id + Rds ‘ ‘ 
Gj +k =1lor2, = 1, 2). 
k ) 


a 
i 
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If we put 


and if we suppose that 


#1 


k l (jk) plik) (kj) (jk) (jk) 
a; = 5 la Loy), = 5 la — V- 


(ik) l (jk) Wk) (kj) jk) / 9 

(jk) (jk) ne 
where b;)"’ andd;, are different from zero when; # k, but equal to zero when 7 


= k, then the solutions (26) become 


t i (10 (20) ‘ 
= fall” cos vt + sin + alt + aff? cos 2 


+ sin 2v = 1, 2), 


“il 


t 0) 0) . Zyl 20) ‘ 
=e“ cosvt + sinvt] +e [5 + ci” cos 2v t 
+ di sin 2v 
The succeeding steps of the integration are similar to the preceding step, and 


an induction to the general term will now be made to show that the process of 


integration may be carried on for any desired number of steps. 


Let us suppose that u,,, vj, 7 = 1, 2, have been computed for / = 0, ..., 
n — 1, and that 
(jk) + 
= > a) e 
(28) i,k 0 
G+k=1,2,...,141). 
0 


where and af’; BY” and BY”? are conjugate complex if #k but real if j=k. 
In order to make the induction it will be necessary to show from the differential 
equations in #;, and v,, obtained from (S’) that the solutions for these variable 
are the same as (28) if / HN. 

Step n + 1: Coefficients of ein (S’). Let us consider the differential equations 


obtained by equating the coefficients of ¢" in (S’) after the various solutions in 


(28) have been substituted. ‘They are found to have the same form as (22) if the 


(27) 
win 
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second subscript on the variables and the right members is changed ton. The 


right members, however, have the form 


U A (7k) ~~ (re 4 
‘lin ’ 
(29) jk=0O0 


i,k 0 


where AY” and BY" are complex quantities similar to those in (28). The com- 
plementary functions of these differential equations which have the desired ex- 
ponentials are the same as (16) if the second subscript on the variables and the 
superscript on the constants of integration are replaced by 7. The particular 
integrals of these differential equations can be found as at the previous step by 


the method of the variation of parameters. They are 


n+ 1 AY + k) n+ 1 
= in + (jk) + 
i,k=0 Al + kysz) | j,k =O 
n+l (Gi +k) n+ 1 
Asn —(jdo + kds)t (jk) + 
= e = Bin e ’ 
Fae Al — + | ik=0 
(30) AY + k) + 
an (ide + (jk) (jr2 + kds)t 
ji k=O A| + kysz) | j,k =0 
+1 (Gj +k) n+ 
Ain (jdz + >» e + kds)t 
"2 ) 2n 
j,k =0 Al — + kds) | j,k =0 


G+tk = 2,3,...,2+ 1). 


The various A’s in the preceding equations are defined as follows: 


A[ — (je + kds)] denotes the determinant A in (10) if D is re 
placed by — (jdA2 + kd;). Since 7 + k is not less than 2, these de- 
nominators do not become A( — de) or A( — ds) and are therefore different 


from zero. 

A,” *) i = 1, 2, 3, 4, denote the preceding determinant if the elements of 
the ith row, reading from top to bottom, are replaced by AY, BY”, AY” and 
BY”, respectively. 


On determining the constants of integration a,” and a,”’ by the initial condi- 


tions (1S), and then unifying the notation as at step 2, the complete solutions 


— 
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which have the desired form are found to be the same as (28) if] = . ‘This 
completes the induction. 


The solutions, therefore, at the general step 7 are 


n+ 
(jk) (jr 
Min Qin ’ 
j,k 0 
> n +1 
Cin Pin ¢ 
i,k 0 


= 1,2. x= O, 1, 2, ). 


When equations (31) are substituted in (15), the asymptotic solutions of the 


differential equations (S) are then found to be 


(jk) + on 
€ it; Qin €, 
€v; = > Bi € 3,2; +8 = I, 2, 
n=1jk=0 


(ik) (jk) 
where the second subscript on the constants a; and ;)' has been made to 
conform with the powers of e. 
If these solutions are expressed in trigonometric form by the substitution 


used in (27) we obtain 


2 2 
Zk ) é 
L(+ ¥ 3, t) > ) * [iain cos(2k — 1) vt 


1 k=Oorl 
333) n « 
2 1, 2h ~~ 1) yal 22) Ip 
+ bi sin(2k — 1) ¢ (Qj, cos 2k vit 
+ sin 2k v 
g:(+ 73, — = 1, 2), 


where S’ denotes that the highest value of 2] —1 or 2/ is n, and =* denotes that 
the lowest values of 2k — 1 and 2k are | and 0, respectively. The function 
g,( + V3, — t) isthe same as f; (+ 1/3, — t) if the constants a and b with their 
subscripts and superscripts are replaced by c and d, respectively, with the same 
subscripts and superscripts. When these equations are substituted in (7) we 
obtain the parametric equations of the orbits which approach the vertices of the 


equilateral triangles in configuration I, viz., 


| 
nn+i 
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l 


/3 
= ~ (te + ms) + + V3, — 2), 


te => my + ms + fe( + V3, — 0) 


(34) 
/3 
m, + ge( + V3, — 0), 
(1, X1 + Me 
M3 
l 
w= (11, Vi + Mle Ve). 
Ms 


The orhjts which approach the configuration IT are obtained by changing the 
sign of V3 in (34). They have the equations 


| 
m=; (m3; — Me) + fi (— 3, — t), 


3 
(m2 + m3) + V/'3, — J), 


m, + m3 + fe(— V3, — 0, 


(35) 
V3 

= + (— — 2), 

| 
(m1 + Me 

My 

Vs = — — + Me 


Kquations (34) and (35) therefore represent the orbits which approach the 
vertices of the equilateral triangles in configurations I and II, respectively, as 
the time approaches plus infinity. These orbits deal with the future of the system 
and, at the risk of being censored, we shall cite them as future I and future U1, 
respectively. 

(. The asymptotic solutions in ¢ 
asymptotic orbits which approach the vertices of the equilateral triangles as 
the time approaches — ”. Such orbits deal with the past while the previous 


orbits forecast the future. So then the mathematical astronomer can say with 


‘and ¢”. ‘This section deals with the 


more than poetic license 


“Backward, turn backward, O time in your flight.” 


= 
=] 
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It is obvious that these orbits can be constructed by making the same use of 
the exponentials e and ec” as was made of e~™ and e~™ in constructing 
the previous orbits. This construction will not be considered in detail as we 
shall show that the past orbits can be obtained from the future orbits in a very 
simple way. ‘This method, it will be observed, is the converse of the adage of 
“history repeating itself.” 

Let us consider the differential equations (8) and let their solutions (33) be 
denoted by f;,( + — 1), + 73, — = 1,2. Next consider the effect 
of changing the signs of t, 4/3, v7; in (8). ‘The left members of the first and third 
equations of (8) remain unchanged while those of the second and fourth change 
signs. Since <_" and o) are even in 1/3, v1, and v, considered together, the above 
changes of signs will leave these expressions unaltered. But as the right members 
of the second and fourth equations in (8) contain the factor V3, these changes of 
signs will produce a change of sign not only in the left members but also in the 
right members of these two equations and the minus sign can be cancelled off 
in both equations. Hence the differential equations (8) are unchanged if the 
signs of t, 1/3, v: and v2 are changed. ‘These changes of signs have no effect upon 
the initial conditions (18) and consequently if we make the same changes of signs 
in the solutions f,( + 4/3, — t), ¢)( + +/3, — 20), it will still leave them solutions 
of (8). Thus f,( — 3, + t) and — ¢g,( — W3, + 2) are solutions of (8) and 
the corresponding solutions of (1) are 


(ms — m) + fi(— V3, + 0), 


"n= v3 (m2 + ms) — g1(— V/3, + 2), 


ty 


ts 
II 


; m + ms + fo(— V3, + 0), 
(36) 
m, — go(— V3, + 2), 


l 

= — — (m + me 
Me 

= — — (mn + me 
Ms 


‘These are the orbits which approach the vertices of the equilateral triangle in 


configuration I as the time approaches — “, and will be cited as past I. The 
corresponding orbits for configuration II, past II, are obtained from (36) by chang- 
ing the sign of +/3. They are 
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x, = (ms Me) + (+ V/3, + t), 


V/3 
— + m) (+ V3, 4+ 9, 


II 


X= 5m + ms + fo(+ V3, + 2), 
(37) 


3 
+ m, — go(+ V3, + 2), 
(m, x1 + me 
M3 
] 
= — — + me 


From the form of the equations of these four orbits it is evident that the past 
II orbits are obtained from the future I orbits by changing the sign of ¢ in the 
latter and reflecting in the x-axis. The same relation exists between past I and 
future IT orbits. 

7. The convergence of solutions. Only the formal construction of the solu- 
tions (8) has been made, and we shall now consider their convergence. 

The convergence of the solutions (32) of the differential equations (8) de- 
pends upon the form of the characteristic exponents in terms of which the so- 
lutions are expanded. Now it has been shown by Poincaré* that such solutions 
as (32), (34) and (35) will converge as t approaches -++ # provided that the real 
parts of the characteristic exponents in terms of which the solutions have been 
expanded are different from zero and positive, and likewise solutions such as 
(36) and (37) will converge as t approaches — provided that the corresponding 
exponents are different from zero and negative. The characteristic exponents 
of this problem are 0, 0, = Ay, = Ae, + Az, equations (13), but only A» and As 
have their real parts different from zero, and it is in terms of only these exponents 
that the solutions have been expanded. Hence the solutions (34) and (35) 
converge as ¢ approaches + @, and the solutions (36) and (37) converge as f 
approaches — 

8. Illustrative examples. We shall conclude this paper with illustrative 
numerical examples. The values of the masses chosen are not the ratios of 
the masses of the Sun, Jupiter and the planetoids mentioned in § 1, but more sim- 
ple values to illustrate the nature of the orbits. 

Let m, = 0.2, me = 0.3 and mz; = 0.5. Then M, the sum of the masses, is 


unity, and #, the mean angular motion, is + | or— 1. We have chosen the 


* Loc. cil., vol. 1, pp. 338-343. 
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counter-clockwise direction of rotation and put » = + 1. ‘The vertices of 
the equilateral triangle for configuration I are (0.1, 0.44/3), (0.6, — 0.14/3) 
and ( — 0.4, — 0.14/3), and for configuration II (0.1, — 0.44/3), (0.6, 0.14/3) 
and ( — 0.4, 0.14/3) for 11), m2 and ms, respectively. 

The computation has been carried out for only the linear terms in ¢ for the 


»). The values of the various constants together with the equations 


orbits (52). 
which define them are listed in the following table. ‘The conjugates of the var 


ious terms are omitted. 


TABLE I 
Constant Equation Value 
| 
A (11 D2(D§+2D1+3 .0925D? +2 .0925 
(12 
7 (14) et seq 0.152+0.948 /—1 
(14) et seq —().246+0.547 
(14) et seq —().630—0.315 
a, a (16), (19 0. 500— (0.225+0.91l4y) 
— Be” (16), (19 — (0.289+0.8667) —(0.440—0.1397) /—1 
a, (16), (19 0.5007 +(0.328+0.225y7) 
a (16), (19 (0.386 +0. 288+) +(0.015—0.578y 


The solutions for € 2), € 7, € #2, and € 7 are 


| cos 0.987t — (0.450 + 1.8287) sin O.987t] + 
e— 0.685 [(0.578 + 1.752y) cos 0.9871 + (0.880 — 0.278y) sin O.987t] 4 
e—0.688¢ | cos + (0.656 + $50y) sin O.987t | 
€ e— 9.6886 [(0.772 + cos + (0.080 — 1.152y) sin O.987t] + 


Example I. 
If we put « = y = U.1 and consider only the linear terms in e¢, these solu- 


tions become 


e—9.688¢ 1000 cos O0.987t — 0.0633 sin 
€i — ¢—9.688/ [0 0751 cos 0.987t + 0.0852 sin 0. 
e—9.6881 [0100 cos 0.987t + sin 0.9877), 
e— 9.688! [0.0830 cos 0 — 0.0085 sin 0.9871}. 


The values of the above quantities for various values of ¢ are listed in Table II. 


rhe corresponding values for ¢ u; and € v7; have been computed from the center 


| 
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of gravity equations. The coérdinates (e im, € 


note the x- and y-displacements of #72), #12 and mis, 


of the equilateral triangle in configuration I. 


Ui), Le 


OS3 
O76 
O70 
063 
O56 
O49 
043 
O36 
O30 
O25 
O19 
O10 
003 
003 
OO7 
O10 
013 
O12 
O09 
O06 
0033 
QOOS 
0014 
00138 


TABLE II 

€ 4 0.1 
0 i} +0.100 | —O.075 || +0.010 +0) 
0.1 +()_ O87 —0. 078 || +0.016 
0.2 +0 O75 —0 O79 i +0. O21 | +0 
0.3 +0063 079 +() O24 | 
0.4 +0052 —0.078 || +0.028 +0 
0.5 +() O41 — O76 } +0.030 +-() 
0.6 +0 032 —(0 +0 031 +0) 
0.7 +0023 069 +().032 +() 
+0. 015 —(0.065 +0033 +() 
0.9 +0.008 | —0.061 +0.033 | +0 
+0.001 | —O.057 || +0.0382 +O 
1.2 |} —0.009 | —0.047 || +0.030 +0) 
1.4 —O.017 | —0.037 || +0.027 +() 
1.6 |} —0.021 —0.028 || +0.023 —() 
1.8 —0.024 | —0.020 |} +0.019 | —O0 
2 —0.025 —(.012 || +0.015 | —O 
2.4 —.022 —(0 001 +0. 008 | —O 
2.8 —O.017 | +0.006 +0). O02 —() 
2.2 —(0 O11 008 —(0.001 —( 
3.6 —0.006 +() —().003 
| —0.0015 | +0.0072 —0.0037 | —O 
1.5 +0.0015 +0 OO46 —0.0032 | 
5 | +0. 0027 +0.0021 —0.0021 +() 
5.5 +0.0026 +0. 0008 —(0.0011 +0 
6 +0). 0019 +0. 0006 00038 +0 
7 +0 +0. 0009 +0.0004 +() 


0005 


Mo, € Ve) and uz, cde 


—(0 046 
044 
—(0.042 
—() 040 
037 
034 
—0 O32 
—( O29 
— O26 
—0 O23 
—0 020 
—(0 O14 
—(0 010 
—0 005 
—0 O02 
+0. O01 
+0. O04 
+0) OOS 
+0005 
+0. 004 
+0 
+0) OO11 
+0. 
—(0 


respectively, from the vertices 


—0 
O15 
—0 O10 
—0 O06 
—() OOD 
+-() OO) 
+0) 
+0 OOS 
+ O10 
+() 
+O O18 
+O O13 
O13 
O12 
+) O11 
+() OOS 
+-() OOD 
+0 OU2 
+0 
—0 OO15 
—0 OO14 
— 0.00097 
—0.00052 
+-() OOOOS 


The diagram of the orbit of m, for Example I is shown in Figure 3. This 


orbit is, of course, with respect to the rotating axes. 


directions of motion. 


Example IT. 


We have also computed the orbits when « = 


0.1 and y = 1. 


The arrow indicates the 


The equations 


determining the displacements, in so far as the linear terms in € are concerned, are 


= e— 9.6881 (0.100 cos 0.9871 
= — e—9.688! (().231 cos 0.9871 
= e— 0.6881 (0.100 cos 0.9871 
= e— 9.688! (0.135 cos 0.9871 


+ 


0.228 
0.060 
0.111 

0.112 


sin 0.987%), 
sin 0.9871), 
sin 0.9871), 


sin 0.9871). 


On substituting in these equations the various values of ¢ as in Table II and mak- 


ing use of the center of gravity equations to determine ¢€ u; and € 73, the following 


displacements are found: 


120 
| 
a 
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* 4 
\ 
Ms \ m2 
Fic. 3 
| \ 
| \ 
\ 
A \) 
1 Mz 


Fic. 4. 
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TABLE III 
e=0.1 v =1 

0 +0.100 —0.231 +0.100 +0.135 || —0.100 | +0011 
0.1 +0 .072 —0.221 +0.103 | +0.115 |} —0.091 | +0.019 
0.2 +0 .046 —0.208 +0.105 | +0.096 || —0.081 +0 .026 
0.3 +0.024 | —0.195 +0.104 | +0.079 || —0.072 | +0.031 
0.4 +0.004 —(0.179 +0.103 | +0.062 | —0.063 +0.034 
0.5 —0.014 —0.164 +0.100 | +0.047 —0.054 | +0.038 
0.6 —0.029 —0.149 +0.096 +0.033 || —0.046 | +0.040 
0.7 —0.042 —0.133 || +0.092 +0.020 —0.038 +0.041 
0.8 —0.053 —0.119 +0 .086 +0.008 —0.031 +0 .042 
0.9 —0.061 —0.104 +0.080 —0.001 —0.024 +0.042 
1.0 —0.068 —0.090 +0.074 —0.010 |} —0.018 +0.041 
1.2 —0.076 | —0.063 +0.062 —0.023 —0.007 +0.039 
1.4 —0.078 —0.039 +0.049 —0.032 +0.002 +0.035 
1.6 —0.076 —0.019 +0.037 —0.038 +0.009 +0.030 
13 —0.071 —0.003 +0 .026 —0.040 +0.013 +0.025 
2 —0.063 +0.009 +0.016 —0.040 +0.016 +0.020 
2.4 —0.044 +0.024 || +0.001 —0.034 +0.017 | +0.011 
2.8 —0.026 +0.028 —0.008 —0.024 +0.015 +0.003 
33 —0.011 +0.026 —0.011 | —0.015 +0.011 —0.001 
3.6 —0.0001 | +0.0198 |} —0.0114 | —0.0066 || +0.0069 —0.0039 
4 +0.0061 | +0.0130 |} —0.0095 | —0.0008 |} +0.0033 | —0.0047 
4.5 +0.0087 | +0.0054 || —0.0061 | +0.0032 || +0.0001 | —0.0041 
5 | +0.0078 | +0.0003 || —0.0028 | +0.0045 || —0.0015 | —0.0028 
5.5 +0.0054 | —0.0024 || —0.0004 | +0.0039 || —0.0019 —0.0014 
6 +0.0028 | —0.0031 |} +0.0009 | +0.0027 || —0.0017 | —0.0004 
7 | —0.0004 | —0.0018 || +0.0012 | +0.0004 || —0.0005 —0.0005 


The diagram of the orbits of the three bodies in this example is found in Fig- 


ure 4. 


length. 


The length of the side of the equilateral triangle is taken as the unit 


In conclusion, the author wishes to express his thanks to his former colleague, 
F. M. Wood, M.A., B.Sc., for verifying the algebraic and numerical computa- 


tions in this paper. 
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Page 51. J. F. Rirt, Prime and composite polynominals. 


Page 63. Line 7, for y-'(v) read y—(w); 
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